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Ozone  (Oa)  is  a  key  species  in  the  troposphere  for  several 

reasons.  It  is  a  key  ingredient  of  "Los  kngeles-type  smog”  as  a 

'  / 

result  of  photocheaical  reactions.  Tropospheric  Oa  also  inpacts  the 
radiation  budget  and  nay  cause  climate  change  because  it  is  a  green- 

l>  * 

house  gas  with  a  strong  absorption  band  centered  at  9.6  /pm  in  the 

atmospheric  window  region  (Fishman  et  al.f  1979).  Also,  photolysis  of 

r ! 

Oa  by  radiation  of  wavelengths  less  than  318  nm  is  the  primary  source 
of  hydroxyl  (OH)  radicals  in  the  tropospheres 
Oa  +  hv  — >  0(1D)  +  Oa  (1) 

0(‘D)  +  M  — >  0{3P)  +  M  (2) 

0 ( 1 D)  +  HaO  -->  2  • OH  (3) 

OOD)  is  an  excited  state  of  atomic  oxygen;  0(aP)  is  the  ground  state. 
OH  radicals  are  important  in  the  cleansing  of  the  atmosphere  through 

oxidation  of  pollutant  species  such  as  carbon  monoxide  (CO) ,  methane 

(CH« ) ,  non-methane  hydrocarbons  (HMHCs)  and  some  halocarbons  (Levy, 
1971). 


Sources  of  tropospheric  Oa  have  been  under  much  debate  in  the 
past  few  years.  Oa  is  known  to  be  photochemically  active  in  the 
stratosphere  where  it  strongly  absorbs  ultraviolet  radiation,  causing 


the  reversal  in  lapse  rate  there.  As  explained  by  Fishman  (1985), 
earliest  studies  considered  Oa  inert  ttn  the  troposphere;  deposition  at 
the  surface  was  believed  balanced  My  injection  from  the  stratosphere 
(Junge,  1962).  Photochemical  production  of  Oa  in  the  troposphere  was 
first  considered  after  the  discovery  of  Oa  as  a  key  species  in  "Los 


Angeles-type  smog”.  In  the  early  1970's,  the  idea  that  significant 
amounts  of  O3  may  be  produced  outside  the  polluted  planetary  boundary 
layer  (PBL)  was  proposed  (Chameides  and  Walker,  1973 ? ^  Crutzen,  1973). 
It  is  now  quite  clear  that  photochemistry  and  transport  both  effect 
the  observed  O3  distribution  in  the  troposphere.  What  remains  is  to 
determine  the  relative  contributions  of  each  process  under  specific 
conditions. 


Though  there  are  various  sources  of  tropospheric  03 ,  by  comparing 
correlations  between  Os  and  other  trace  species,  and  considering  the 
meteorological  situation,  one  can  often  draw  conclusions  about  the 
origin  and  transport  of  O3  in  the  troposphere. 


Carbon  monoxide  has  its  largest  source  in  the  PBL,  from  combus¬ 
tion  of  fossil  fuels  or  biomass  burning.  In  addition,  CO  is  the  end 
product  in  the  oxidation  chain  of  CH«  and  NMHCs.  The  primary  loss 
mechanism  for  CO  in  the  troposphere  is  reaction  with  the  OH  (Levy  II, 
1971;  Condon  et  al.,  1987).  In  mid-latitudes,  the  typical  lifetime  of 


CO  is  of  the  order  of  a  few  months  (Thompson  and  Cicerone,  1986).  CO 
was  found  to  be  a  good  tracer  of  PBL  air  (Fishman  et  al.,  1980; 
Pickering  et  al.,  1988a)  so  when  CO  and  O3  concentrations  are  simul¬ 
taneously  high,  the  source  of  the  elevated  O3  is  the  PBL,  or  in-situ 
production  in  the  free  troposphere  by  PBL  precursors  (discussed 
below).  CO  concentration  decreases  sharply  in  the  stratosphere,  thus 
high  O3  with  low  CO  indicates  a  stratospheric  source  of  Os . 


Like  CO,  water  vapor  (HsO)  originates  in  the  PBL,  at  the  surface, 


jpci.  lal 


ropy 

,*Speor co 


-V  «•  .V  V.  -\ 


v-rvjys 


ORIGIN  OF  TROPOSPHERIC  OZONE 
OVER  CENTRAL  US 

by 

Donna  P.  McNamara 

Thesis  submitted  to  the  Faculty  of  the  Graduate  School 
of  the  University  of  Maryland  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 
Master  of  Science 
1988 

Advisory  Committee: 

Dr.  Russell  Dickerson 
Dr.  Robert  Ellingson 
Dr.  George  Huffman 


ACKNOWLEDGEMENTS 

I  would  like  to  thank  ay  husband,  Brien,  for  his  constant  help 
and  support  throughout  ny  time  in  graduate  school.  Be  also  aided  my 
research  by  calculating  zenith  angles  and  local  times,  and  copying 
photographs . 

1  would  also  like  to  thank  Dr.  Russell  Dickerson,  my  research 
advisor,  for  his  help  selecting  a  topic,  and  his  constant  support  and 
encouragement.  He  showed  great  patience  and  persistence  teaching 
chemistry  to  a  meteorologist,  and  was  always  available  when  I  needed 
help. 


I  also  received  much  support  from  my  fellow  graduate  students  who 
collected  the  field  data.  Ken  Pickering  provided  innumerable 
assistance,  and  was  like  a  second  research  advisor.  Winston  Luke  and 
Linda  Nunnermacker  helped  open  my  eyes  to  analytical  chemistry  and 
explained  the  workings  of  the  various  trace  gas  measurement 
techniques.  Brian  Gockel  prepared  the  framework  for  the  numerical 
model  used  to  estimate  ozone  production;  He  and  Melody  Owens  helped 
me  stay  ahead  in  my  many  battles  with  the  computer. 

Dr.  George  Huffman  provided  very  helpful  support  and  ideas  for 
the  meteorological  analysis  in  this  research. 

Jose  Meitin,  of  NOAA/ERL,  promptly  provided  much  supporting  data, 
such  as  rawinsondes,  radar  images,  and  satellite  pictures. 


ft 


Si 


fl 

i 


a 

.N 


ii 


r; 

f, 

* 


lW*S 


Dr.  Arlan  Krueger,  of  NASA,  provided  the  initial  TOMS  pictures 
which  sparked  my  interest  in  the  TOMS  analysis.  Dr.  Jack  Fishman, 
also  of  NASA,  supplied  the  enhanced  imagery  and  raw  data  used  in  the 
analysis.  The  analysis  was  aided  by  P.K.  Bartia,  of  STX  Corp.,  who 
explained  the  processing  algorithms. 

The  Research  Aviation  Facility  of  the  National  Center  for 
Atmospheric  Research  supported  the  aircraft  measurements. 

This  research  was  supported  by  NSF  Grants  ATM-83-05843  and  ATM- 
86-19491,  and  the  PRECP  and  PRESTORM  programs. 

Finally,  thanks  must  go  to  my  father,  Wallace  Palmer,  for  his 
great  help  with  the  word  processing  and  printing  of  this  thesis. 


TABLE  OF  CONTENTS 


List  o f  Tables 
List  of  Figures 
Introduction 

Chapter  I  Instrumentation 
Chapter  II  Synoptic  Situation 
Chapter  III  Results 

A.  Flight  1 

1.  Vertical  Profiles 

2.  Time  Series 

B.  Flight  2 

1.  Vertical  Profiles 

2.  Time  Series 
Chapter  IV  Discussion 

1.  Correlations  -  Flight  1 

2.  Correlations  -  Flight  2 

3.  Event  1  -  High  O3  and  CO  at  5.5  km 

4.  Event  2  -  Large  Os  Gradient  at  10.6  km 

5.  TOMS  Data 

Chapter  V  Tropospheric  Ozone  Production 
Chapter  VI  Summary 
Appendix  A  -  Data  Procesing 
Appendix  B  -  1  Minute  Average  Data 


References 


LIST  OF  TABLES 


Number 

Page 

•V 

> 

1. 

Overall  Trace  Gas  Correlations 

50 

1  i*. 

2. 

Trace  Gas  Profile  Correlations 

58 

| 

3. 

Ozone  Production  at  Altitudes  of  Level  Flight 

91 

4. 

Reactions  Involved  in  Ozone  Production 

92 

1 

5. 

Hydrocarbon  Analysis 

97 

•V 

.V 

6. 

Effective  Methane  Calculation 

100 

A 

> 

»•*: 

7. 

Net  Ozone  Column  Production 

102 

I 


> 


:>i 

sH 


hi 


i 


u>. 


I 


pJIM IW  W  ICTKSk?  W*  V “  >r.  •HT'CV't;  wvv  yj yj yv wv; \r;y.  v. ■«-»«-.'  v 


LIST  OF  FIGURES 

Number 

1.  300  nb  heights  for  1200  GMT,  22  June  1985. 

2.  Sea  level  isobars  and  frontal  positions  for  1200  GMT 
22  June  1985. 

3.  GOES  IR  imagery,  1200  GMT,  22  June  1985. 

4.  Flight  track  for  the  first  mission. 

5.  Cross  section  of  the  first  flight  track. 


6. 

7. 


GOES  IR  imagery,  2300  GMT,  22  June  1985. 

Flight  track  for  the  second  mission. 

8.  Sabreliner  altitude  versus  time  for  the  second  flight. 

9.  Vertical  profile  of  Oa  and  CO,  ascending  leg,  flight  1. 

10.  Vertical  profile  of  Oa  and  CO,  descending  leg,  flight  1. 

11.  Vertical  profile  of  temperature  and  dee  point,  ascending 
leg,  flight  1. 

12.  Vertical  profile  of  temperature  and  dew  point,  descending 
leg,  flight  1. 

13.  Vertical  profile  of  NO,  NOx,  and  NOy ,  flight  1. 

14.  GOES  IR  imagery,  1801  GMT,  22  June  1985. 

15.  Time  series  of  O3  and  CO,  flight  1. 

16.  Time  series  for  NO,  NOx  and  NOy,  flight  1. 

17.  Time  series  of  dew  point  depression,  flight  1. 

18.  GOES  IR  imagery,  1900  GMT,  22  June  1985. 

19.  GOES  IR  imagery,  2200  GMT,  22  June  1985. 

20.  Vertical  profiles  of  O3  and  CO,  all  data,  flight  2. 

21.  Average  vertical  profiles  of  O3  and  CO,  flight  2. 


Number 


22.  Vertical  profile  of  temperature  and  dew  point,  ascending 
leg,  flight  2. 

23.  Vertical  profile  of  temperature  and  dew  point,  descending 
leg,  flight  2. 

24.  Vertical  profile  of  NO,  NO* ,  and  N0» ,  flight  2. 

25.  Time  series  of  O3  and  CO,  flight  2. 

26.  Time  series  for  NO,  NO*  and  NOy ,  flight  2. 

27.  Time  series  of  dew  point  depression,  flight  2. 

28.  Scatter  plot  of  CO  versus  O3,  all  data,  flight  1. 

29.  Scatter  plot  of  CO  versus  O3  in  the  free  troposphere, 
flight  1. 

30.  Scatter  plot  of  NO*  versus  Os,  all  data,  flight  1. 

31.  Scatter  plot  of  NO  versus  03,  all  data,  flight  1. 

32.  Scatter  plot  of  NOy  versus  O3 ,  all  data,  flight  1. 

33.  Vertical  profiles  of  O3 ,  CO  and  TD,  ascending  leg, 
flight  1. 

34.  Vertical  profiles  of  O3 ,  CO  and  TD,  descending  leg, 
flight  1. 

35.  Scatter  plot  of  CO  versus  O3 ,  all  data,  flight  2. 

36.  Scatter  plot  of  CO  versus  O3  in  levels  II  and  Ha, 

flight  2. 

37.  Scatter  plot  of  CO  versus  Os  in  level  III,  flight  2. 

38.  Scatter  plot  of  CO  versus  O3  in  level  IV,  flight  2. 

39.  Scatter  plot  of  NO,  NOx  and  N0y  versus  O3  at  10.6  km, 
flight  2. 

40.  Scatter  plot  of  NOx  versus  O3  at  10.6  km,  flight  2. 


Nuaber 


41.  Scatter  plot  of  NO  versus  Os,  all  data,  flight  2. 

42.  Average  vertical  profiles  of  O3 ,  CO  and  TD  for  flight  2. 

43.  Flight  track  of  Sabreliner  at  10.6  km  on  flight  2. 

44.  NO*  and  NOy  versus  O3 ,  after  the  first  peak  in  O3  at 
10.6  ka,  flight  2. 

45.  Hind  speed  versus  Os  at  10.6  ka,  flight  2. 

46.  200  mb  heights  for  1200  GMT,  21  June  1985;  with  back 

trajectories  of  air  saapled  at  10.6  km  on  flight  2. 

47.  TOMS  imagery  from  the  Nimbus  7  satellite  for  21  June  1985 

48.  TOMS  imagery  from  the  Nimbus  7  satellite  for  22  June  1985 

49.  Gridded  TOMS  data  from  22  June  1985. 


INTRODUCTION 


Ozone  (Os)  is  a  key  species  in  the  troposphere  for  several 
reasons.  It  is  a  key  ingredient  of  "Los  Angeles-type  smog"  as  a 
result  of  photochemical  reactions.  Tropospheric  O3  also  impacts  the 
radiation  budget  and  may  cause  climate  change  because  it  is  a  green¬ 
house  gas  with  a  strong  absorption  band  centered  at  9.6  jim  in  the 
atmospheric  window  region  (Fishman  et  al.,  1979).  Also,  photolysis  of 
O3  by  radiation  of  wavelengths  less  than  318  nm  is  the  primary  source 
of  hydroxyl  (OH)  radicals  in  the  troposphere: 

O3  +  hv  — >  0(1D)  +  O2  (1) 

OOD)  +  M  -->  0(3P)  +  M  (2) 

0(XD)  +  H2O  — >  2  OH  (3) 

0{1D)  is  an  excited  state  of  atomic  oxygen;  0(3P)  is  the  ground  state. 
OH  radicals  are  important  in  the  cleansing  of  the  atmosphere  through 
oxidation  of  pollutant  species  such  as  carbon  monoxide  (CO) ,  methane 
(CH« ) ,  non-methane  hydrocarbons  (NMHCs)  and  some  halocarbons  (Levy, 
1971). 

Sources  of  tropospheric  O3  have  been  under  much  debate  in  the 
past  few  years.  O3  is  known  to  be  photochemically  active  in  the 
stratosphere  where  it  strongly  aLsorbs  ultraviolet  radiation,  causing 
the  reversal  in  lapse  rate  there.  As  explained  by  Fishman  (1985), 
earliest  studies  considered  O3  inert  in  the  troposphere;  deposition  at 
the  surface  was  believed  balanced  by  injection  from  the  stratosphere 
(Junge,  1962).  Photochemical  production  of  O3  in  the  troposphere  was 
first  considered  after  the  discovery  of  O3  as  a  key  species  in  "Los 


Angeles-type  snog".  In  the  early  1970's,  the  idea  that  significant 
amounts  of  O3  nay  be  produced  outside  the  polluted  planetary  boundary 
layer  (PBL)  was  proposed  (Chaneides  and  Walker,  1973;  Crutzen,  1973). 
It  is  now  quite  clear  that  photochemistry  and  transport  both  effect 
the  observed  Oa  distribution  in  the  troposphere.  What  remains  is  to 
determine  the  relative  contributions  of  each  process  under  specific 
conditions. 

Though  there  are  various  sources  of  tropospheric  Oa ,  by  comparing 
correlations  between  Oa  and  other  trace  species,  and  considering  the 
meteorological  situation,  one  can  often  draw  conclusions  about  the 
origin  and  transport  of  O3  in  the  troposphere. 

Carbon  monoxide  has  its  largest  source  in  the  PBL,  from  combus¬ 
tion  of  fossil  fuels  or  biomass  burning.  In  addition,  CO  is  the  end 
product  in  the  oxidation  chain  of  CH«  and  NMHCs.  The  primary  loss 
mechanism  for  CO  in  the  troposphere  is  reaction  with  the  OH  (Levy  II, 
1971;  Condon  et  al.,  1987).  In  mid-latitudes,  the  typical  lifetime  of 
CO  is  of  the  order  of  a  few  months  (Thompson  and  Cicerone,  1986).  CO 
was  found  to  be  a  good  tracer  of  PBL  air  (Fishman  et  al.,  1980; 
Pickering  et  al.,  1988a)  so  when  CO  and  O3  concentrations  are  simul¬ 
taneously  high,  the  source  of  the  elevated  O3  is  the  PBL,  or  in-situ 
production  in  the  free  troposphere  by  PBL  precursors  (discussed 
below).  CO  concentration  decreases  sharply  in  the  stratosphere,  thus 
high  Oa  with  low  CO  indicates  a  stratospheric  source  of  Oa . 

Like  CO,  water  vapor  (HaO)  originates  in  the  PBL,  at  *-he  surface, 


but  it  is  not  an  ideal  tracer  because  precipitation  is  a  sink  for  B2O. 


The  sources  of  nitrogen  oxides  are  more  complex.  The  largest 
source  of  nitric  oxide  (NO)  is  combustion  of  fossil  fuels  and  biomass 
burning.  But  another  large  source  is  lightning,  and  smaller  amounts 
of  NO  come  from  soil  exhalation,  stratospheric  intrusions,  aircraft 
emissions,  oxidation  of  ammonia  (NHa)  and  photolysis  of  marine  N02~ 
(Stedman  and  Shetter,  1983;  Logan,  1983).  Estimates  of  the  relative 
contributions  of  each  source  vary  dramatically.  For  example,  Davis  et 
al.  (1987)  estimate  lightning  to  contribute  from  2  to  20  Tg  N/yr. 
Because  the  sources  are  numerous  and  variable,  and  only  recently  have 
effective  measurement  techniques  been  developed,  the  estimates  of 
typical  NO  concentrations  for  specific  situations  are  also  highly 
uncertain.  It  is  common  to  consider  NO  and  nitrogen  dioxide  (NO2) 
together  as  NOx  and  all  odd  nitrogen  (NO,  NO2 ,  NO3 ,  N2O0,  HNO2 ,  HNOa , 
HO2NO2,  peroxyacetyl  nitrate  (FAN)  and  other  organic  nitrates)  as  NOy . 
In  the  troposphere  NOx  and  O3  participate,  in  part,  in  the  following 
reactions: 

NO  +  03  — >  NO2  +  O2  (4) 

NO2  +  hv  — >  NO  +  0  (5) 

0  +  O2  +  M  —  >  03  +  M.  (6) 

These  reactions  approach  a  steady-state,  called  the  photostationary 
state,  when  not  disturbed  by  changes  in  light  intensity,  reactant 
concentrations,  or  competing  reactions  (Leighton,  1961).  When  NO2  is 
photolyzed  (by  radiation  of  wavelengths  less  than  420  nm) ,  O3  will 
eventually  be  produced.  If  there  is  enough  NOx  in  the  presence  of  CO, 
CH< ,  NMHCs  and  sunlight,  the  peroxy  radicals  (HO2  and  R02)  formed  will 


convert  NO  to  NOs  without  consuming  Oa  (or  0) ,  resulting  in  an 
increase  in  O3  (Chameides  and  Walker,  1973;  Crutzen,  1973).  Thus, 
production  of  03  is  critically  dependent  on  NOx  concentration  (e.g. 
Crutzen  and  Gidel,  1983;  Liu  et  al.,  1987).  The  concentration  of  NOx 
nay  vary  by  as  much  as  4  orders  of  magnitude  from  the  remote  marine 
atmosphere  to  the  polluted  continental  PBL  (Logan,  1983).  NOx  has  a 
short  lifetime,  just  one  to  two  days  in  summer,  before  it  is  converted 
to  one  of  the  NO?  species,  such  as  nitric  acid  (HNO3),  which  can  act 
as  a  reservoir  or  be  removed  from  the  atmosphere.  NOy  is  lost  through 
wet  deposition  of  soluble  species  such  as  HNO3 ,  and  dry  deposition  of 
NOs  and  HNOs  (Logan,  1983). 

Tropospheric  concentrations  of  a  trace  gas  are  in  part  a  function 
of  its  lifetime  and  sources,  as  well  as  the  meteorological  conditions. 
Concentrations  of  trace  gases  such  as  oxides  of  nitrogen  exhibit  large 
variations  in  the  troposphere  because  they  react  on  the  time  scale  of 
minutes  to  hours,  their  sources  are  not  uniformly  distributed,  and 
they  participate  in  chemical  reactions  before  meteorological  transport 
can  uniformly  distribute  them  throughout  the  troposphere.  Chemistry 
as  well  as  meteorology  must  be  considered  when  analyzing  the  distribu¬ 
tion  of  these  species.  In  contrast,  trace  gases  such  as  CH4  have  very 
long  lifetimes,  so  transport  has  time  to  mix  these  species  thoroughly, 
resulting  in  concentrations  that  are  nearly  constant  throughout  the 
troposphere.  Intermediate  to  these  are  species  such  as  CO,  which  do 
show  variability  in  concentration  in  the  troposphere,  but  have 
lifetimes  long  enough  to  be  considered  unreactive  on  a  time  scale  of  a 
few  hours;  location  of  sources  and  transport  are  then  responsible  for 


most  of  the  observed  variations. 


Vhen  pollutants  remain  in  the  PBL,  they  are  generally  short  lived 
due  to  wet  and  dry  deposition.  Wind  speeds  are  lower  in  the  PBL,  so 
transport  is  minimized.  If  the  pollutants  are  transported  to  the  free 
troposphere  however,  their  lifetimes  increase,  and  they  can  be 
transported  much  farther  (Dickerson  et  al.,  1987).  Pollution  now 
becomes  a  regional  or  national  problem,  rather  than  a  local,  primarily 
urban,  problem.  Recent  research  by  Liu  et  al.  (1987)  suggests  the 
efficiency  of  Os  production  increases  as  NOx  concentration  decreases; 
potentially  more  Oa  is  produced  than  current  estimates  if  more  NO* 
escapes  from  the  PBL  and  is  diluted  in  the  free  troposphere.  Gener¬ 
ally  the  PBL  is  topped  with  a  temperature  inversion  that  inhibits 
mixing  with  the  free  troposphere,  and  tropospheric  layers  persist 
because  mean  vertical  velocities  on  the  synoptic  scale  are  of  the 
order  of  just  a  few  centimeters  per  second  (Holton,  1979).  However, 
when  instability  allows  cumulus  convection,  mesoscale  vertical  mixing 
may  be  large.  Gidel  (1983)  and  Chatfield  and  Crutzen  (1984)  developed 
photochemical  models  that  included  transport  by  cumulus  clouds;  their 
results  suggest  convective  transport  from  the  PBL  plays  a  larger  role 
in  tropospheric  chemistry  than  previously  believed.  Isaac  et  al. 
(1984)  found  evidence  that  precipitating  and  non-precipitating  clouds 
transport  pollutants  out  of  the  PBL.  Chameides  et  al.  (1987)  found 
model-calculated  OH  concentrations  were  correlated  with  measured 
elevated  NO  concentrations  from  a  thunderstorm  source  and  concluded 
that,  as  the  convection  pumps  PBL  air  to  high  altitudes,  the  anvil 


area  of  a  thunderstorm  should  be  a  "photochemical  hotspot". 


During  the  summer  of  1985,  tropospheric  trace  gases  and  meteoro¬ 
logical  data  were  collected  as  part  of  a  project  to  study  trace  gas 
profiles  over  rural  North  America.  The  project,  "Trace  Gas  Profiles 
II”  was  sponsored  by  the  National  Science  Foundation  (NSF) ,  and  was 
conducted  concurrently  with  two  other  projects,  PRE-STORM  (sponsored 
by  the  National  Oceanic  and  Atmospheric  Administration  (NOAA) ,  and 
NSF)  and  PRECP  (sponsored  by  the  US  Dept,  of  Energy).  Operations  were 
based  in  Oklahoma  City.  The  National  Center  for  Atmospheric  Research 
(NCAR)  Sabreliner  aircraft  was  fitted  with  detectors  to  measure  the 
trace  gases  ozone,  carbon  monoxide,  and  nitrogen  oxides;  as  well  as 
many  meteorological  parameters.  This  thesis  examines  two  Sabreliner 
flights  on  June  22,  1985.  For  analysis  of  other  flights  from  this 
project  see  Dickerson  et  al.  (1987),  Pickering  (1987),  Pickering  et 
al.  (1988a),  and  Pickering  et  al.  (1988b). 

Dickerson  et  al.  (1987)  studied  a  thunderstorm  and  found  elevated 
values  of  CO,  NMHCs,  O3 ,  NOx  and  NOy  in  the  outflow  region  of  the 
storm,  suggesting  thunderstorms  are  a  mechanism  for  significant 
vertical  transport  from  the  PBL.  In  contrast,  Pickering  et  al. 
(1988a) ,  found  a  thunderstorm  case  where  a  cold  front  cut  off  the  cell 
from  PBL  air,  such  that  the  cell  could  not  mix  PBL  air  with  the  free 
troposphere.  Therefore,  significant  vertical  transport  cannot  be 
assumed  with  all  convective  cells. 

The  two  flights  considered  here  investigated  convection  associ¬ 
ated  with  a  dissipating  frontal  system  and  fair  weather  in  the  region 


u1 


later  the  same  day.  The  purpose  of  this  research  is  to  present  the 
observed  trace  gas  distributions,  account  for  the  origin  of  gradients 
in  tropospheric  O3 ,  and  quantify  the  Os  production  possible  under  the 
conditions  observed.  Vertical  profiles  and  time  series  of  trace  gases 
are  presented.  Correlations  between  the  various  trace  gases  are 
analyzed.  Total  Ozone  Happing  Spectrometer  (TOMS)  data  from  the  NASA 
Nimbus  7  polar-orbiting  satellite  are  presented  to  look  for  effects  of 
tropospheric  O3  concentrations  on  column  abundances.  An  estimate  of 
the  in-situ  photochemical  production  of  O3  is  made  using  a  simple 
photochemical  model. 


I.  INSTRUMENTATION 


Data  were  collected  froa  the  NCAR  Sabreliner-- a  twin  engine 
research  jet  with  a  flight  altitude  ceiling  of  about  11  km.  In 
addition  to  its  normal  package  of  meteorological  equipment  (NCAR, 
1984),  the  Sabreliner  carried  instruments  to  measure  Oa ,  CO,  NO,  NOx 
and  NOy .  Details  of  the  trace  gas  data  processing  are  described  in 
Appendix  A.  Detailed  one  minute  average  data  for  both  flights  are 
given  in  Appendix  B.  For  some  missions,  grab  samples  of  NMHCs  were 
collected  for  later  analysis. 

Nitrogen  oxides  were  measured  by  O3  chemiluminescence  with  a 
Thermo  Electron  Corporation  Model  14B  detector  modified  for  greater 
sensitivity  (Delany  et  al.,  1982;  Dickerson  et  al.,  1984).  Measure¬ 
ment  of  NO  using  chemiluminescence  is  based  on  reaction  4,  above,  of 
the  photostationary  state,  in  which  NO  and  O3  react  to  form  N02 
(Stedman  et  al.  1972).  Approximately  10  percent  of  the  N02  formed 
will  be  in  an  excited  state  and  will  emit  light,  which  can  be  mea¬ 
sured.  Before  oxides  of  nitrogen  other  than  NO  can  be  detected,  they 
must  first  be  converted  to  NO.  The  detector  uses  ferrous  sulfate  to 
convert  NOx  to  NO  and  a  molybdenum  catalyst  to  convert  NOy .  An 
equipment  intercomparison  by  Fehsenfeld  et  al.  (1987)  determined  there 
is  a  positive  interference  in  the  measured  NOx  due  to  PAN  when  using 
the  ferrous  sulfate  converter.  To  obtain  measurements,  the  sample 
flow  is  mixed  with  O3 .  A  background  signal  is  obtained  by  sending 
sample  air  through  a  preconverter  where  the  NO  and  O3  react  out  of 
view  of  the  photomultiplier  tube.  The  detector  mode  was  usually 


changed  every  3  minutes  to  measure  NO,  NO*,  NOy  and  background  on  a 
12-minute  cycle.  A  measurement  was  recorded  every  second  and  averaged 
over  10  second  intervals;  about  30  seconds  of  data  were  lost  during 
each  mode  change.  The  instrument  has  a  detection  limit  of  about  20 
ppt  NO;  the  uncertainty  of  NO  detection  is  +  14  percent,  and  that  of 
NOy  is  t  20  percent  at  the  95  percent  confidence  level  (Fehsenfeld  et 
al.  1987). 

Ozone  was  measured  using  a  Dasibi  Model  1003  detector,  which 
measures  O3  absorption  of  254  nm  radiation.  A  reference  signal  is 
determined  by  passing  the  sample  air  first  through  a  catalyst  which 
destroys  all  the  Os .  The  sample  and  reference  signals  are  subtracted 
to  obtain  the  O3  concentration  every  7  seconds.  The  instrument  has  a 
detection  limit  and  uncertainty  of  about  +  2  ppb  (at  95  percent 
confidence  level) . 

Carbon  monoxide  was  measured  using  a  modified  Thermo  Electron 
Corporation  Model  48  gas  correlation  non-dispersive  infrared  detector. 
Gas  filter  correlation  (GFC)  analyzers  operate  by  sensing  the  absorp¬ 
tion  from  CO  in  the  4.67  jin  fundamental  vibration-rotation  band. 
Radiation  in  this  region  is  also  absorbed  by  the  wings  of  other  trace 
gases,  most  significantly  carbon  dioxide  (COs)  and  HsO.  GFC  corrects 
for  this  by  passing  the  source  beam  alternately  through  a  filter  of 
pure  nitrogen  (Ns)  and  a  filter  of  pure  CO,  before  passing  through  the 
sample  cell.  The  COs  and  HsO  in  the  sample  will  attenuate  both  beams 
almost  the  same  amount,  but  CO  in  the  sample  cell  will  change  only  the 
beam  that  passed  through  the  Ns .  The  CO  cell  produces  a  reference 


signal.  The  difference  in  the  intensities  of  the  two  beams  yields  the 
concentration  of  CO  in  the  sample  cell.  Modifications  further 
eliminated  the  HsO  interference  by  removing  it  from  the  sample  with  a 
drying  tube,  reduced  the  amount  of  zero  drift  by  inserting  a  chemical 
zero,  and  increased  the  signal-to-noise  ratio.  The  modified  instru¬ 
ment  has  a  detection  limit  of  about  24  ppb  (S5  percent  confidence 
level)  with  a  60  second  signal  average  (Dickerson  and  Delany,  1988). 
Values  were  recorded  every  10  seconds,  and  the  instrument  was  usually 
in  zero  mode  for  3  of  every  12  minutes.  Data  used  in  this  research 
are  30  second  signal  averages. 


II.  Synoptic  Situation 


On  21  June  1985  an  upper  level  trough  was  located  in  the  mountain 
states,  extending  as  far  south  as  southern  Colorado  at  300mb. 
Associated  with  this  trough  was  a  strong  cold  front  located  north  of 
Oklahoma,  through  central  Nebraska  and  along  the  northern  border  of 
Colorado.  Around  0200  GMT  (2100  CDT)  convection  developed  along  a 
squall  line  east  of  Amarillo  TX.  By  0500  GMT  a  squall  line  also 
developed  along  the  border  of  Kansas  and  Oklahoma.  Throughout  the 
night,  the  line  pushed  to  the  southeast.  In  Kansas,  there  were  three 
reports  of  tornadoes;  hail  and  severe  wind  gusts  were  reported  in 
Oklahoma  and  Kansas  (Meitin  and  Cunning,  1985) . 

By  the  morning  of  22  June,  the  upper  level  trough  (Figure  1)  had 
become  vertically  stacked  through  the  central  US  from  Minnesota  to 
Missouri,  and  continued  to  recede  northward  in  the  next  24  hours.  The 
surface  front  at  1200  GMT  (Figure  2)  extended  from  Illinois  to  Texas, 
through  central  Oklahoma,  and  had  weakened  drastically.  GOES  IR  (NOAA 
geostationary  infrared)  satellite  imagery  (Figure  3)  showed  high  level 
clouds  remaining  from  central  Oklahoma  to  southern  Illinois  along  the 
old  front.  Radar  reports  confirmed  some  convection  remained, 
organized  in  a  line  that  pushed  slowly  through  central  Arkansas  during 
the  flight  periods.  The  first  flight  of  the  day  took  off  from 
Oklahoma  City  at  1220  CDT  (1720  GMT)  and  flew  a  nearly  straight  east- 
southeast  out  and  back  track  through  central  Arkansas,  turning  around 
over  northwest  Mississippi  (Figure  4) .  The  Sabreliner  ascended 
quickly  to  an  altitude  greater  than  8  km,  and  spent  much  of  the  flight 
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in  cloud  (Figure  5) .  The  return  leg  wee  in  clear  air  after  the 
descent  to  5.5  km  just  before  reaching  the  Arkansas-Oklahoma  border. 
The  flight  sampled  a  vertically  mixed  troposphere,  shortly  after  the 
peak  in  convective  activity  as  determined  by  satellite  imagery. 

By  2300  GMT  (1800  CDT) ,  the  remaining  convection  continued  to 
decay  (Figure  6) .  The  surface  front  had  dissipated  just  north  of 
Oklahoma  City.  The  second  flight  of  the  day  left  Oklahoma  City  at 
1702  CDT  and  flew  east,  making  several  loops  within  Oklahoma  (Figure 
7).  The  Sabreliner  spent  much  time  in  level  flight  at  7  km  and  10.6 
km,  and  made  an  extra  penetration  into  the  PBL  (Figure  8) .  The  only 
clouds  were  fair  weather  cumulus,  built  mostly  to  about  2  km,  with 
some  as  high  as  4.5  km. 

All  day,  the  area  from  Oklahoma  City  through  eastern  Arkansas 
stayed  predominantly  under  southerly  flow  at  the  surface  from  the  Gulf 
of  Mexico,  never  experiencing  an  obvious  frontal  passage.  At  upper 
levels  (above  700mb) ,  the  wind  had  a  westerly  component  throughout  the 
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Fig.  6.  GOES  IR  satellite  imagery,  2300  GMT  (1800  CDT) , 
22  June  1985. 
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Ill .  RESULTS 


A.  flight  1  (1215-1455  CDT) 


Analysis  of  the  trace  gas  data  fron  the  first  flight  shows  strong 
evidence  of  vertical  transport  of  pollution  out  of  the  PBL  over  large 
areas,  possible  in-situ  photochemical  production  of  ozone  and  strong 
vertical  mixing  throughout  the  troposphere. 


1.  Vertical  Profiles 

Vertical  profiles  of  Os  and  CO  concentration  for  the  ascending 
(profile  1)  and  descending  (profile  2)  legs  are  presented  in  Figures  9 
and  10.  The  profiles  for  flights  1  and  2  show  high  levels  of  O3  and 
CO,  and  substantial  variability  when  compared  to  remote  clean  air 
profiles  (e.g.  Danielson  et  al.,  1987;  Fishman  et  al.,  1987)  measured 
over  the  Pacific  Ocean.  Temperature  and  dew  point  are  also  plotted 
against  altitude  in  Figures  11  and  12. 


The  intense  solar  radiation  of  early  afternoon  in  summer  produced 
conditions  favorable  for  convection,  making  it  hard  to  determine  the 
exact  height  of  the  PBL.  For  this  research,  the  PBL  will  be  assumed 
to  be  capped  by  the  slight  temperature  inversion  at  about  1.9  km  on 
both  ascent  and  descent,  which  best  matches  the  first  drop  in  CO  that 
appears  in  profile  2  (this  definition  of  PBL  includes  layers  that 
could  also  be  called  the  mixed  and  cumulus  layers).  Maximum  CO 
concentration  in  the  PBL  was  almost  160  ppb,  which  is  relatively  clean 
compared  to  the  average  nonurban  surface  value  for  that  latitude  of 
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about  200  ppb  reported  by  Logan  et  al.  1981.  O3  concentration  was 
fairly  low  in  the  PBL,  ranging  from  about  20  ppb  at  take  off  to  40  ppb 
at  the  end  of  the  flight. 

A  surprising  feature  is  the  strong  increase  in  both  O3  and  CO 
above  the  PBL  up  to  about  5  km;  this  feature  was  also  evident  during 
the  late  afternoon  flight.  Concentration  of  both  species  began  to 
decrease  right  above  the  temperature  inversion  at  4.7  km  in  profile  1 
and  4.5  km  in  profile  2.  The  CO  concentrations  in  this  layer  ap¬ 
proached  PBL  values,  03  concentrations  were  much  higher  than  in  the 
PBL,  and  the  air  was  dry.  Unfortunately,  no  NO  data  are  available  for 
this  layer.  Figure  13  shows  odd  nitrogen  concentrations  versus 
altitude  measured  for  the  entire  flight.  NOx  values  up  to  1400  ppt 
were  recorded  in  the  PBL,  but  a  secondary  maximum  (approximately  600 
ppt)  was  evident  just  below  4  km.  This  indicates  more  than  enough  N0X 
to  catalyze  the  photochemical  production  of  O3  (Crutzen  and  Gidel, 
1983,  estimate  that  greater  than  5  ppt  NO  at  the  surface  and  20  ppt 
aloft  will  result  in  O3  production) . 

The  Oa  concentrations,  higher  than  those  in  the  PBL,  may  have 
resulted  from  a  combination  of  in-situ  production  and  venting  of  PBL 
air.  With  a  PBL  1.9  km  thick,  even  with  a  dry  deposition  velocity  of 
1  cm/s  (Droppo,  1985) ,  the  O3  lifetime  in  the  PBL  would  be  more  than  2 
days;  O3  could  escape  the  PBL  before  being  lost  at  the  surface. 
Alternatively,  this  layer  could  be  the  original  PBL  that  was  lifted 


and  displaced  by  convection-induced  downdrafts  of  cleaner  mid-tropo- 
spheric  air  (e.g.,  Lyons  et  al.,  1986a).  It  is  important  to  note  that 


one  can  determine  air  of  PBL  origin  (by  its  elevated  CO,  KOx ,  NKHC  or 
HaO) ,  but  as  the  PBL  was  not  sampled  over  large  areas  beneath  the 
flight  track  and  upstream,  one  can  make  no  conclusions  about  the 
original  geographic  location  of  the  PBL  air.  The  air  may  have 
traveled  long  horizontal  distances. 

Above  this  layer,  the  vertical  profiles  indicate  widely  varying 
concentrations  of  O3 ,  CO  and  nitrogen  oxides.  There  is  some  layering 
present,  but  it  is  not  consistent  between  the  ascending  and  descending 
profiles.  The  large  variations  in  the  horizontal  probably  resulted 
from  the  strong  convective  mixing  on  this  day.  These  results  are 
consistent  with  those  of  Lyons  et  al.  (1986)  and  Lyons  and  Calby 
(1986) ,  who  describe  large  vertical  transport  by  mesoscale  systems 
ranging  in  size  from  single  airmass  thunderstorms  to  large  (up  to  2500 
km  wide)  convective  complexes.  The  Sabreliner  was  often  in  cloud  at 
these  higher  altitudes  (Figure  5).  Analysis  of  radar  reports, 
satellite  imagery  (Figure  14) ,  and  nose  and  side  camera  films  indicate 
some  regions  of  active  convection  and  some  regions  of  outflow  cirrus. 
At  1700  GMT,  National  Weather  Service  (NWS)  radar  at  Little  Rock, 
Arkansas  indicated  a  line  of  showers,  about  50  km  wide,  extending 
east-northeast  0  west-southwest  through  central  Arkansas.  The 
showers  were  dissipating,  but  some  remained  as  the  Sabreliner  crossed 
the  line  around  1800  GMT.  Winds  from  the  west  aloft  extended  the 
cirrus  outflow  to  the  east.  In  general,  above  5  km  CO  ranged  from 
about  70  to  140  ppb,  O3  was  between  30  and  70  ppb,  and  odd  nitrogen 
showed  large  variation. 


Fig.  14.  GOES  IR  satellite  imagery,  1801  GMT  (1301  CDT) , 
22  June  1985. 


2.  Tine  Series 

The  Sabreliner  spent  a  significant  amount  of  time  at  constant 
altitude,  so  time  series  were  prepared  for  the  trace  gas  species. 
Figure  15  shows  O3  and  CO  versus  time;  Figures  16  and  17  show  odd 
nitrogen  and  dew  point  depression  versus  time  (dew  point  depression  is 
the  difference  between  the  dew  point  and  ambient  temperature;  it  shows 
good  variation  at  all  altitudes,  unlike  other  moisture  parameters  such 
as  mixing  ratio  which  drop  to  near  zero  out  of  the  PBL) .  Only 
constant  altitude  data  are  included  in  these  figures. 

A  very  interesting  feature  shows  up  in  the  time  series  at  7000 
sec  (1418  CDT) ,  lasting  for  a  little  more  than  6  minutes  (about  60  km 
of  flight  trajectory) .  The  CO  climbed  to  150  ppb  and  the  O3  to  its 
highest  value  of  the  flight,  93  ppb.  The  Sabreliner  was  in  level 
flight  at  5.5  km  in  clear  air  returning  to  Oklahoma  City.  The  air  was 
very  dry,  as  indicated  by  the  large  dew  point  depression.  Satellite 
imagery  at  1800  GMT  and  1900  GMT  (Figures  14  and  18)  show  the  convec¬ 
tion  was  rapidly  deteriorating  in  Arkansas  and  western  Oklahoma  before 
this  time.  The  nose  camera  film  showed  nothing  visibly  different 
about  the  air  as  the  concentrations  increased;  uniform  cumulus  cells 
topping  below  flight  level  were  scattered  in  the  foreground.  Unfor¬ 
tunately,  the  nitrogen  oxide  detector  was  being  calibrated  at  this 
time.  The  origin  of  this  air  is  discussed  later  as  Event  1. 

Another  surprising  feature  is  seen  in  the  time  series  for  the  odd 
nitrogen  data.  The  NO  concentrations  generally  run  well  below  500 
ppt,  except  at  2500  sec.  and  3200  sec.  (about  1804  and  1816  GMT)  where 
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Fig.  18.  GOES  IR  satellite  imagery,  1900  GMT  (1400  CDT) 
22  June  1985. 


they  reach  760  and  960  ppt.  The  Sabreliner  was  in  cloud  at  both  times 
at  flight  levels  above  10  km.  The  NO*  and  NOy  values  were  not 
unusually  high.  The  Sabreliner  may  have  encountered  NO  from  a  recent 
lightning  strike  or  a  fresh  aircraft  plume;  NO  would  be  high,  but  may 
not  have  yet  been  substantially  converted  to  NO*  or  the  other  NOy 
species.  Radar  screen  images  (plan  position  indicator,  250  nm)  from 
the  NWS  at  Little  Rock  at  1802  GMT  confirm  the  existence  of  a  line  of 
thunderstorms  over  Little  Rock  extending  east-northeast  through  west- 
southwest,  moving  southward.  Winds  aloft  carried  the  cirrus  shield 
toward  the  east.  It  is  impossible  to  quantify  exactly  the  distance 
from  a  cell  because  the  Sabreliner  was  almost  directly  over  Little 
Rock  at  1804  GMT,  and  the  antenna  was  set  at  an  elevation  angle  which 
passed  the  radar  beam  under  the  convection  close  to  the  radar  site. 
Nevertheless,  lightning  is  possible  within  the  anvil  itself,  even 
after  the  thunderstorm  cell  has  decayed  (US  Air  Force,  1982). 
Chameides  et  al.  (1987)  also  attributed  elevated  NO  concentrations 
above  500  ppt  to  lightning  while  near  a  storm  anvil. 

B.  Flight  2  (17:02-19:12  CDT) 

The  mission  flown  in  the  afternoon,  under  fair  skies,  is  an 
excellent  contrast  to  the  case  described  above.  Analysis  of  trace  gas 
data  showed  distinct  layers  immediately  above  the  PBL  and  in  the  upper 
troposphere,  with  pollutant  concentrations  nearly  as  high  as  in  the 


1.  Vertical  Profiles 


Satellite  Imagery  near  takeoff  time  (Figure  19)  shows  do  remain¬ 
ing  organized  cumulus  convection  in  eastern  Oklahoma.  Vertical 
profiles  for  Os  and  CO,  which  include  all  data,  are  shown  in  Figure 
20.  Except  for  small  fluctuations,  the  Os  profiles  overlay  each  other 
quite  well.  Where  CO  data  exist  more  than  once  at  an  altitude,  the 
values  are  also  quite  similar.  The  atmosphere  appears  quite  homogen¬ 
eous  in  the  horizontal  and  layered  in  the  vertical.  Therefore,  it  is 
reasonable  to  consider  O3  and  CO  profiles  obtained  by  averaging  all 
available  data  to  0.1  km  layers  (Figure  21).  The  temperature  and  dew 
point  profiles  (Figures  22  and  23)  helped  delineate  the  layers  found 
in  Figure  21.  Odd  nitrogen  profiles  are  shown  in  Figure  24. 

In  the  PBL  (I)  below  2.1  km,  CO  and  NOx  decrease  with  height  and 
Os  increases  with  height.  CO  ranges  from  about  130  to  175  ppb,  NOx 
from  about  500  to  1500  ppt,  and  O3  from  about  20  to  50  ppb.  Above  the 
PBL  up  to  about  4.1  km  is  once  again  a  very  dirty  layer  (II)  charac¬ 
terized  by  drier  conditions,  and  O3 ,  CO  and  NOx  concentrations 
reaching  as  high  as  75  ppb,  172ppb,  and  1200  ppt,  respectively.  Nose 
and  side  camera  footage  show  scattered  fair  weather  cumulus  mostly  to 
about  2  km,  but  some  penetrating  the  lowest  temperature  inversion  and 
reaching  about  4.5  km.  The  dirty  air  in  layer  II  may  be  the  top  of  an 
old  PBL  that  was  well-mixed  to  the  higher  level  and  cut  off  from  the 
surface  when  a  new  inversion  formed  near  the  surface  at  night  (Vuko- 
vich  et  al.,  1985).  Concentrations  in  layer  II  could  then  increase 
from  in-situ  photochemical  O3  production  and  further  venting  of  PBL 


Fig.  19.  GOES  IR  satellite  imagery,  2200  GMT  (1700  CDT) 
22  June  1985. 
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Fig.  20.  Vertical  profiles  of  03  and  CO,  all  data,  fligh 
2.  03  is  7  sec.  data.  CO  is  30  sec.  averages,  with  break 
in  data  where  the  instrument  was  in  background  mode 
Constant  altitude  data  are  not  included. 
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flight  2.  All  available  data  were  averaged  in 


TEMP  and  DEW  POINT  vs.  ALTITUDE,  PROFILE  1, FLIGHT  2 


Stull  (1984)  parameterized  cumulus  clouds  as  forced  clouds  and 
active  clouds.  Forced  clouds  resemble  the  shallow  clouds  described 
above,  which  only  mix  air  within  the  PBL  and  provide  conditions  for 
aqueous-phase  chemical  reactions.  Active  clouds,  such  as  the  taller 
clouds  mentioned  above  which  penetrate  the  temperature  inversion 
capping  the  PBL,  vent  PBL  air  to  higher  altitudes  and  increase 
pollutant  levels  aloft.  Layer  II  is  also  topped  by  a  temperature 
inversion  beginning  at  about  4.2  km,  which  prevents  further  upward 
mixing  and  allows  pollutant  levels  to  accumulate.  Pollutant  levels 
decrease  rapidly  above  the  inversion  up  to  4.6  km  (layer  Ha) . 

Layer  III,  between  4.6  and  7.8  km,  contains  the  cleanest  air 
sampled  during  the  mission.  The  layer  was  dry  and  none  of  the  fair 
weather  cumulus  were  seen  penetrating  to  this  layer.  CO  concentra¬ 
tions  were  as  low  as  65  ppb,  Oa  concentrations  decreased  to  28  ppb 
(the  only  lower  concentration  was  found  near  the  surface),  and  NOx 
concentrations  were  mostly  less  than  500  ppt.  Above  7.8  km  (layer 
IV) ,  pollutant  levels  were  variable  but  generally  higher  than  in  layer 
III.  Oa  concentrations  exceeded  70  ppb  (as  high  as  141  ppb  in  level 
flight  at  10.6  km),  CO  concentrations  reached  120  ppb,  and  NOx 
concentrations  were  variable.  Back  (constant  level)  trajectory 
analysis  of  the  air  at  9.1  km  (30,000  ft)  placed  this  air  in  north 
Texas  and  northeast  New  Mexico  when  convection  was  occurring  in  that 
region.  As  summertime  convection  to  the  upper  troposphere  is  common, 
and  probably  carries  Oa  or  Oa  precursors  to  the  upper  troposphere,  we 
might  expect  this  layer  of  the  atmosphere  to  show  elevated  concen- 


trations  of  CO  and  NOx . 


2.  Time  Series 

Time  series  were  prepared  for  data  from  this  flight  as  well, 
again  including  data  only  from  level  flight  segments.  Figure  25  shows 
CO  and  Oa,  Figure  26  shows  odd  nitrogen,  and  Figure  27  shows  dew  point 
depression  versus  time. 

The  most  striking  feature  of  all  the  time  series  is  in  the  O3 
time  series  at  10.6  km  (beginning  at  3400  sec.).  The  first  jump  in 
O3 ,  from  less  than  40  ppb  to  more  than  110  ppb,  was  such  a  surprise 
that  the  onboard  scientist  requested  the  pilot  fly  a  loop  and  reinves¬ 
tigate  the  same  location  again.  The  existence  of  the  feature  was 
confirmed,  and  Os  concentration  this  time  reached  more  than  140  ppb. 
Although  at  other  levels  many  smaller  peaks  in  Os  match  peaks  in  CO 
(1500  sec.,  2300  sec.,  2900  sec.,  and  6200  sec.),  these  two  peaks  at 
10.6  km  are  independent  of  changes  in  CO.  The  origin  of  this  sharp 
gradient  is  discussed  later  as  Event  2. 
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5,  time  series  of  dew  point  depress 
egments  of  level  flight,  flight  2. 


IV.  DISCUSSION 


1.  Correlations  -  Flight  1 

To  understand  further  the  relative  importance  of  tropospheric  03 
production  versus  stratospheric  injection,  one  must  make  a  more 
detailed  assessment  of  the  correlations  between  the  various  trace  gas 
species  (Fishman  et  al.  1987).  One  way  to  do  this  is  to  prepare 
scatter  diagrams  between  the  individual  measurements  of  two  species, 
but  this  technique  has  problems.  The  average  or  background  profiles 
of  the  different  species  vary  significantly  and  will  overpower  the 
small-scale  correlations  between  the  species.  For  example,  water 
vapor  tends  to  be  very  large  in  the  PBL  and  becomes  very  small  above; 
CO  is  also  largest  in  the  PBL,  but  shows  a  much  more  gradual  decrease 
aloft;  Oa  tends  to  be  large  in  the  upper  troposphere  due  to  its 
stratospheric  source;  and  odd  nitrogen  tends  to  show  two  maxima,  in 
the  PBL  and  upper  troposphere.  In  addition,  it  is  very  hard  to  prove 
a  statistical  correlation  between  the  two  species  because  each 
measurement  in  a  time  series  is  not  independent  and  the  degrees  of 
freedom  are  much  less  than  the  total  number  of  measurements. 

The  effective  degrees  of  freedom  (N')  for  a  time  series  can  be 
computed  according  to  the  following  formula: 

N'  ■  N  /  (1  +  rin'  +  r*r2'  +  ^rs’  +  r«r4'  +  roro') 

N  is  the  number  of  data  pairs;  ri ,  r* ,  etc.  are  autocorrelation 
coefficients  for  the  first  variable  in  the  series;  n',r2',  etc.  are 
autocorrelation  coefficients  for  the  second  variable  in  the  series 
(Quenouille,  1952).  Prior  to  computing  the  correlations,  the  data 


were  paired  and  times  were  matched  by  taking  the  species  with  the 
longest  measurement  interval  and  averaging  all  measurements  of  the 
second  species  occurring  more  than  half  in  the  longer  time  interval 
(for  example,  a  CO  30  second  average  could  match  up  with  as  many  as 
five  Oa  7  second  averages) .  The  correlations  using  ail  data  from  the 
flight  are  listed  in  Table  1.  Most  of  the  regression  coefficients  are 
small,  and  all  are  statistically  insignificant  using  the  effective 
degrees  of  freedom.  However,  scatter  plots  of  some  subsets  suggest 
relationships  between  the  trace  gases. 

CO  versus  O3  (Figure  28)  appear  positively  correlated  with  a 
slope  of  0.24  (units  ppb/ppb) ,  but  this  relationship  is  not  quite 
significant  at  the  95  percent  confidence  level  (R2  of  0.199).  The 
small  number  of  points  with  high  CO  and  Os  suggest  there  may  be  at 
least  places  where  the  two  species  are  well  correlated.  Including 
only  data  above  the  PBL  (Figure  29)  eliminates  the  known  level  of  poor 
correlation  right  at  the  surface  (where  O3  has  its  shortest  lifetime 
due  to  removal  by  dry  deposition  and  chemical  reactions) .  This 
analysis  produces  a  slope  of  0.32,  and  R2  of  0.312  meets  the  signifi¬ 
cance  test  at  the  95  percent  confidence  level  (Panofsky  and  Brier, 
1968) . 

Most  scatter  plots  between  O3  and  the  various  nitrogen  oxides  do 
not  demonstrate  clear  relationships  because  odd  nitrogen  concentra¬ 
tions  show  large  variations;  the  more  interesting  plots  are  shown 
here.  Figure  30  shows  NOx  versus  03;  linear  regression  produces  only 
a  weak,  non-significant  slope  of  -0.0081  (units  ppb/ppt,  R2  of  0.130), 


TABLE  1 


OVERALL  TRACE  GAS  CORRELATIONS 


Flight  1 

Slope 

R* 

Significant 

slope 

CO/Oa 

0.24 

0.20 

no 

CO/DPdep.* 

0.095 

0.080 

no 

03/DPdep. 

0.60 

0.0074 

no 

NO/Oa 

0.0060 

0.062 

no 

NOx /03 

-0.0081 

0.13 

no 

NOy/Oa 

0.0044 

0.12 

no 

Flight  2 

CO/Oa 

-0.17 

0.050 

no 

CO/DPdep. 

-0.041 

0.029 

no 

Oa/DPdep. 

-0.34 

0.016 

no 

NO/Oa 

-0.0028 

0.00028 

no 

NOx /O3 

-0.012 

0.17 

no 

NOy/Oa 

-0.016 

0.026 

no 

*DPdep  -  dew  point  depression 

All  correlations  were  calculated  using  all  data  from  the  flight. 
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seme  regions  of  the  graph  show  positive  correlation  and  some  appear  to 
show  slightly  negative  or  no  correlation.  Some  of  the  positively 
correlated  points  are  between  about  3  and  4.2  km  in  the  inferred  upper 
boundary  layer,  and  the  negatively  correlated  points  include  points 
from  the  PBL  and  upper  troposphere.  Figure  31,  a  plot  of  NO  versus 
Oa,  shows  no  correlation,  which  is  not  surprising  as  the  points  to  the 
far  right  include  the  NO  spikes  possibly  attributable  to  lightning  and 
the  points  with  low  NO  are  near  the  instrument  detection  limit.  The 
plot  of  NOy  versus  Os  (Figure  32)  includes  only  data  above  7.9  km;  the 
slight  positive  slope  of  0.0044  (still  insignificant  with  R3  just 
0.122)  is  expected  because  both  have  sources  in  the  polluted  PBL  air 
and  the  stratosphere.  These  results  differ  from  those  reported  by 
Pickering  (1987),  for  another  fair  weather  case  from  the  same  project, 
where  NO*  and  O3  were  always  positively  correlated. 

Fishman  and  Seiler  (1983)  devised  a  technique  to  quantitatively 
measure  the  correlation  between  the  vertical  fine  structure  of  two 
species  and  overcome  the  problem  of  their  average  profile  variations 
overshadowing  their  small  scale  variations.  Looking  at  concurrent 
vertical  profiles  of  two  species,  they  find  the  best  linear  fit  line 
for  each  profile  and  subtract  this  line  from  the  data  averaged  for 
every  0.1  km  layer.  They  then  compute  the  product  of  the  residuals  at 
each  level  where  data  exist  for  both  species.  The  species  are 
positively  correlated  in  layers  in  which  the  product  is  positive,  and 
negatively  correlated  where  the  product  is  negative.  The  products  are 
summed  over  the  entire  profile  ana  normalized  by  dividing  by  the  sum 
of  the  products  of  the  absolute  values  of  the  residuals.  If  this 
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Scatter  plot  of  NO  versus  O3 ,  all  data,  flight 
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quotient  is  +  1.0  the  species  are  perfectly  positively  correlated,  if 
-  1.0  they  are  perfectly  negatively  correlated;  a  result  of  zero 
indicates  no  correlation.  Fishman  and  Seiler  (1983)  arbitrarily  chose 
0.5  to  be  the  lower  limit  for  significance,  and  conclude  that  bands  of 
elevated  O3  and  CO,  which  occur  with  elevated  NOx  and  NMHCs,  strongly 
suggest  in-situ  photochemical  production  of  O3  is  occurring. 

Fine  structure  correlations  from  this  data  set,  between  the 
possible  combinations  of  O3 ,  CO  and  dew  point  temperature  (TD)  for 
each  profile,  are  summarized  in  Table  2.  The  profiles  are  averaged 
for  each  0.1  km  layer  and  the  best  fit  lines  are  shown  in  Figures  33 
and  34.  The  net  positive  slope  for  CO  in  profile  1,  although  it  only 
contains  data  above  about  5  km,  suggests  how  strongly  the  recent 
convection  disturbed  the  troposphere.  Correlations  with  TD  must  be 
interpreted  with  caution  because  the  lag  time  of  dew  pointer  on  the 
Sabreliner  varies  significantly  with  altitude,  and  is  different  for 
ascent  and  descent.  The  lyman  alpha  instrument  has  almost  immediate 
response,  but  the  values  of  all  moisture  parameters  except  TD  drop  off 
so  abruptly  above  the  PBL  that  they  cannot  be  reasonably  fit  by  a 
linear  regression.  (An  algorithm  exists  from  NCAR  to  calculate 
absolute  humidity  from  the  lyman  alpha  voltages,  but  the  correction  of 
the  lag  is  not  perfect  as  the  algorithm  uses  absolute  humidity 
measured  with  the  dew  pointer  to  calibrate  the  voltages.)  In  addi¬ 
tion,  in  this  case  there  was  a  sink  for  water  vapor  in  the  recent 
rainfall  (2.2  cm  in  the  24  hours  ending  at  23  June  1200  GMT  was 
reported  at  Little  Rock,  for  example)  which  explains  how  CO  could  be 
poorly  correlated  with  TD  when  both  have  sources  in  the  PBL. 


TABLE  2 

TRACE  GAS  PROFILE  CORRELATIONS 


Flight 

Profile 

Species 

Normalized 
Corr.  Coeff. 

Signifii 

1 

1 

03 /CO 

0.0022 

no 

1 

1 

CO/TD 

0.29 

DO 

1 

1 

03  /TD 

-0.26 

no 

1 

2 

Os /CO 

0.48 

no 

1 

2 

CO/TD 

0.063 

no 

1 

2 

Os/TD 

-0.82 

yes 

2 

All 

03 /CO 

0.90 

yes 

2 

All 

CO/TD 

0.61 

yes 

2 

All 

03 /TD 

0.35 

no 
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Like  Figure  32,  vertical  profiles  of  03 ,  CO  and  TD 
descending  leg,  flight  1. 


f 


CO  and  Oa  are  not  highly  correlated;  their  correlation  m 
profile  1  is  near  zero  and  for  profile  2  is  0.48,  positive  but  not 
significant.  Profile  1  shows  high  CO/low  O3  above  10  km  and  at  5  km; 
high  Os/low  CO  at  6.5  km.  The  high  CO/low  O3  regions  could  be  air 
recently  brought  up  by  convection  from  the  PBL  which  would  have  low  Os 
due  to  surface  deposition,  and  low  photochemical  03  production  due  to 
the  cloudy  conditions.  The  high  O3/I0W  CO  layer  could  be  the  result 
of  some  stratospheric  or  upper  tropospheric  air  mixing  downward  and 
being  diluted.  It  is  important  to  remember  (as  Chameides  et  al.,  1987 
point  out)  high  03/low  CO  does  not  automatically  eliminate  the 
possibility  of  in-situ  O3  production;  the  production  of  O3  is  criti¬ 
cally  dependent  on  NO  concentration,  which  has  sources  different  from 
CO.  Profile  2  has  several  spikes  above  8  km  with  high  CO/low  O3 ,  and 
has  high  03/low  CO  between  5.7  and  7.6  km,  which  lowered  its  positive 
correlation.  Certainly  broad  areas  of  the  atmosphere  were  highly  dis¬ 
turbed  by  mesoscale  convective  mixing  during  this  flight.  This  is  an 
area  in  which  photochemistry  is  almost  certain  to  produce  substantial 
amounts  of  O3 ,  but  the  photochemical  production  has  not  yet  had  time 
show  up  in  the  vertical  profiles.  In  Profile  2,  the  Sabreliner 
encountered  less  disturbed  air  that  moved  in  behind  the  convection,  so 
correlation  between  O3  and  CO  was  higher. 


2.  Correlations  -  Flight  2 

As  with  the  first  flight,  correlations  were  computed  between  the 
trace  gases  to  further  understand  the  various  concentrations  observed 
(see  Table  1).  In  Figure  35,  a  plot  of  CO  versus  O3  (all  data),  data 
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Scatter  plot  of  CO  versus  O3,  all  data,  flight 


from  10.6  ko  account  for  the  points  with  high  03  and  little  variation 
in  CO;  and  most  of  the  data  with  high  CO  and  little  variation  in  03 
cone  froa  the  PBL.  Looking  at  subsets  of  the  data  is  more  informa¬ 
tive.  CO  versus  O3  in  just  layers  II  and  Ha  (Figure  36)  yields  a 
positive  slope  of  0.24  (ppb/ppb) ,  but  is  not  significant  with  R*  of 
0.4631  because  this  data  subset  has  only  2.8  effective  degrees  of 
freedom.  In  layer  III,  (Figure  37)  the  slope  is  only  0.014  (also  not 
significant  with  R*  only  0.000526).  In  layer  IV,  (Figure  38)  without 
the  data  at  10.6  km  (which  is  dominated  by  the  strong  Os  peaks),  Os 
and  CO  show  a  higher  positive  correlation  with  slope  0.49  (though  R2 
of  0.546  is  insignificant  due  to  only  2.8  effective  degrees  of  free¬ 
dom)  ,  suggesting  that  some  enhanced  O3  in  the  upper  troposphere  may 
indeed  be  the  result  of  tropospheric  production. 

A  few  subsets  of  odd  nitrogen  versus  O3  suggest  some  correlation. 
At  10.6  km  (Figure  39),  Os  appears  slightly  positively  correlated  with 
NO,  NOx  and  NOy  when  one  ignores  the  large  spikes  in  NO  and  NOy ; 
indeed,  Os  versus  NOx  versus  Os  at  10.6  km  (Figure  40)  has  a  positive 
slope  of  0.206  (ppb/ppt)  (though  insignificant  with  Rs  0.647  and 
effective  degrees  of  freedom  3.3).  Os  with  all  NO  would  be  correlated 
if  one  ignores  the  strong  NO  spike  (Figure  41). 

When  profile  correlations  were  done  using  the  technique  of 
Fishman  and  Seiler  (1983)  discussed  above,  the  results  were  much  more 
straightforward  than  those  for  the  first  flight.  The  air  sampled 
above  about  700  mb  on  this  flight  had  not  been  involved  in  convection 
for  12  hours  or  more,  and  strata  of  Os  had  time  to  develop.  Average 
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Fig.  39.  Scatter  plot  of  NO,  N0X  and  N0y  versus  03  at  10. 
km,  flight  2. 
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profiles  of  O3,  CO,  and  TD  are  shown  in  Figure  42  with  the  best  fit 
line  for  each,  and  the  correlations  are  shown  in  Table  2.  The  TD 
profile  is  a  bit  jagged  above  5  km  as  the  descending  profile  was  drier 
than  the  ascending  one  and  some  0.1  km  layer  averages  contain  more 
points  from  the  ascent  or  descent.  However,  the  levels  of  saturation 
were  the  same  on  ascent  and  descent.  O3  and  CO  were  strongly  correla¬ 
ted  with  a  profile  correlation  coefficient  of  0.90.  Without  precipi¬ 
tation  to  act  as  a  sink  for  water  vapor,  CO  and  TD  were  also  posi¬ 
tively  correlated  (0.61).  03  and  TD  were  also  positively  correlated 
(0.35),  though  not  at  a  significant  level.  This  strongly  suggests 
that  much  of  the  enhanced  03  is  of  tropospheric  origin,  and  that  it 
takes  hours  or  days  after  O3  precursors  are  released  to  the  free 
troposphere  to  produce  correlations  with  0s . 


3.  Event  1  -  High  Os  and  CO  at  5.5  km  (flight  1) 

When  the  Sabreliner  was  returning  to  Oklahoma  City,  it  encoun¬ 
tered  an  area  of  dry  air  with  elevated  Os  (93ppb)  and  CO  (150  ppb) 
concentrations  about  60  km  wide  (see  the  sharp  peaks  at  5.5  km  in 
Figures  15  and  17) .  Convection  in  Arkansas  and  western  Oklahoma  had 
been  rapidly  deteriorating,  and  only  scattered  cumulus  cells  topping 
below  flight  level  remained.  Air  with  this  much  CO  must  have  come 
from  the  PBL,  and  the  dryness  indicates  the  air  must  have  recently 
subsided  nearly  3  km  (determined  simply  from  the  dry  adiabatic  lapse 
rate  of  about  10  degrees  per  km  and  a  dew  point  depression  of  almost 
30  degrees).  Potential  temperature  increased  only  about  a  degree  in 
this  air  over  the  surrounding  values;  equivalent  potential  temperature 
(which  is  a  function  of  the  parcels  temperature  and  moisture,  and 
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accounts  for  the  latent  heat  added  if  all  moisture  were  condensed  out) 
decreased  about  4°C  in  the  dirty  air.  Vertical  wind  velocity  measured 
by  the  Sabreliner  showed  that  on  the  average  the  air  was  sinking  at 
0.15  m/s  (with  a  large  standard  deviation  of  +  0.28  m/s).  The  Oa 
values  are  not  inconsistent  with  air  from  the  upper  troposphere,  but 
alone  give  little  information  on  the  air's  origin  because  the  elevated 
values  could  have  come  from  PBL  air  transported  to  this  altitude  by 
convection,  from  a  stratospheric  intrusion,  or  from  in-situ  production 
by  ozone  precursors.  The  strong  correlation  between  O3  and  CO  indi¬ 
cates  the  Os  could  not  have  come  from  the  stratosphere  (Fishman  et 
al.,  1980).  Because  of  the  recent  convection,  isentropic  back  trajec¬ 
tories  are  not  rel:  ^.le  in  this  case. 

A  likely  scenario  is  the  air  containing  elevated  CO  was  brought 
out  of  the  PBL  by  a  thunderstorm  earlier  that  day,  and  during  the 
rapid  decay  of  the  convection  some  mesoscale  effect  forced  air  to 
rapidly  sink  in  that  small  area.  0a  was  likely  increasing  in  this  air 
parcel  because  PBL  air  contains  the  necessary  precursors  for  in-situ 
photochemical  production.  The  concentration  gradients  would  slowly 
disappear  due  to  eddy  diffusion  with  the  surrounding  air,  but  the 
temperature  gradients  would  disappear  much  more  rapidly,  by  a  mecha¬ 
nism  such  as  gravity  waves,  until  the  air  was  in  thermal  equilibrium 
with  the  surroundings.  This  suggests  that  chemical  tracers>  of 
atmospheric  motion  remain  discernible  longer  than  the  more  traditional 
tracer — potential  temperature. 


4.  Event  2  -  Large  O3  Gradient  at  10.6  km  (flight  2) 

The  Sabreliner  twice  encountered  a  sharp  gradient  in  Oa  while  in 
level  flight  at  10.6  km  (Figure  25).  Figure  43  shows  the  Sabreliner 
flight  track  at  10.6  km.  Maximum  Oa  was  found  to  the  north  and 
minimum  O3  was  found  to  the  south. 

Immediately,  a  PBL  source  or  photochemical  production  source  for 
the  Oa  gradient  had  to  be  eliminated  because  the  O3  concentrations 
were  completely  uncorrelated  with  CO,  which  showed  only  relatively 
small  fluctuations  (about  +  10  ppb)  at  10.6  km.  There  was  no  discern¬ 
ible  change  in  moisture  (HsO  was  below  the  detection  limit  of  the  dew 
pointer  and  lyman  alpha  sensor) .  There  was  a  decreasing  trend  in  the 
odd  nitrogen  data  as  it  cycled  from  NOy  to  NOx  and  NO  as  the  O3 
decreased  after  each  peak.  Plots  of  the  NO*  and  NOy  versus  Os  after 
the  first  peak  are  shown  in  Figure  44.  Positive  correlation  between 
Os  and  odd  nitrogen  with  no  correlation  with  CO  point  to  a  strato¬ 
spheric  source,  but  why  the  sharp  gradient? 

Potential  temperature  did  not  show  much  change  at  10.6  km  (less 
than  2°  C) ,  and  was  uncorrelated  with  the  O3  change,  so  there  was  no 
evidence  the  Sabreliner  had  found  the  edge  of  a  tropopause  fold  from 
an  abrupt  intrusion  of  stratospheric  air.  The  synoptic  situation 
precluded  a  tropopause  fold  because  the  upper  level  trough  was  too  far 
to  the  northeast,  and  in  a  fold  one  would  expect  anticorrelation 
between  Os  and  CO  (Danielson  et  al.,  1987;  Hipskind  et  al.,  1987), 
which  was  not  found  in  this  case. 


40 


35.20N  96.20W 


Fig.  43.  Flight  track  of  Sabreliner  at  10.6  km  on  flight 
2,  with  contours  of  O3  concentration  drawn  from  7  sec. 
data . 


The  only  measured  variable  that  had  a  time  series  matching  the 
change  in  O3  was  wind.  Winds  were  higher  with  a  slightly  more 
southerly,  component  where  03  was  lowest  (Figure  45).  In  fact,  the 
correlation  between  O3  concentration  and  wind  speed  is  significant  at 
the  99  percent  level,  using  the  effective  degrees  of  freedom  for  a 
time  series  as  discussed  above.  An  isentropic  back  trajectory, 
initialized  using  the  Oklahoma  City  rawinsonde  winds,  showed  the  air 
at  10.6  km  had  come  from  the  west  in  the  vicinity  of  the  upper  level 
cut  off  low.  As  suggested  by  the  Oa-wind  correlation,  the  trajec¬ 
tories  were  repeated  for  two  parcels  (Figure  46) ,  one  from  the 
northern  (high  O3)  part  of  the  track  and  one  from  the  southern  (low 
03)  part  of  track,  initialized  with  the  actual  winds  measured  by  the 
Sabreliner  for  each  parcel  for  6  hours  (255  degrees  at  11  m/s  for  the 
northern  parcel  and  245  degrees  at  15  m/s  for  the  southern  parcel) . 
The  two  parcels  quickly  diverged.  The  trajectories  showed  very  little 
change  in  altitude  and  no  interaction  with  convection.  The  air  with 
higher  O3  concentrations  had  come  north  of  the  low  from  near  the  jet, 
while  the  air  with  low  O3  had  come  from  around  the  southern  side  of 
the  low  where  the  tropopause  was  much  higher.  It  appears  the  Sabre¬ 
liner  had  flown  through  a  confluence  zone  between  air  originating  near 
the  jet,  where  the  parcel  had  a  greater  chance  of  encountering  high  O3 
concentrations,  and  air  originating  several  kilometers  below  the 
tropopause  in  a  more  tropical  location,  where  03  concentrations  would 
be  expected  to  be  lower. 

The  isentropic  trajectories  were  constructed  by  hand  using  all 
available  upper  air  data.  Before  the  last  6  hours,  winds  for  the 


Wind  speed  versus  03  at  10.6  km,  flight 


CDT) ,  21  June  1985.  Isentropic  back  trajectories  of  air  sampled  at 
10.6  km  on  flight  2  are  shown,  initialized  with  aircraft  winds  from  the 
north  and  south  of  the  flight  track  for  six  hours.  Box  designates 
geographic  extent  of  griddod  TOMS  data  discussed  in  a  later  section. 


appropriate  potential  temperature  surface  were  determined  for  each 
upper  air  station  using  linear  interpolation  between  winds  observed  at 


200  and  250  mb.  Careful  streamline  analyses  were  done  for  each 
synoptic  observation  time,  and  the  parcel  was  moved  back  along  a 
streamline  using  the  average  horizontal  wind  from  the  surrounding 
stations.  For  example,  streamlines  at  1200  GKT  were  used  to  moved  the 
parcel  from  0600  to  1800  GMT. 

5.  TOMS  Data 

The  TOMS  sensor  onboard  the  NASA  Nimbus  7  satellite  provides 
total  column  O3  (TC03)  over  the  entire  earth  once  a  day.  Fishman  et 
al.  (1986)  have  shown  that  in  the  tropics  many  of  the  features  in  TC03 
are  from  tropospheric  O3  variability,  even  though  tropospheric  O3 
accounts  for  only  5  to  15  percent  of  the  column.  Figure  47  shows  data 
from  June  21,  and  Figure  48  shows  data  from  June  22  (flight  day);  the 
time  of  the  data  is  close  to  local  noon  at  any  location  (the  equator¬ 
ial  crossing  time  is  11:50  am,  in  an  ascending  orbit).  These  data 
were  analyzed  to  determine  if  the  tropospheric  features  found  by 
aircraft,  in  particular  the  high  O3  sampled  at  5.5  km  near  the 
Oklahoma-Arkansas  border  on  the  first  flight  or  the  sharp  O3  gradient 
at  10.6  km  on  the  second  flight,  were  detectable  in  the  TC03  values. 
The  data  are  in  Dobson  units  (Du),  where  one  Du  equals  0.01  mm  of  0» 
if  all  the  O3  in  the  atmosphere  were  compressed  to  standard  tempera¬ 
ture  and  pressure.  The  thickness  corresponding  to  a  typical  TCos 
amount  of  300  Du  is  0.3  cm. 

A  preliminary  look  at  Figure  48  was  encouraging.  The  TCos 
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maximum  (15  Du  above  the  surrounding  area)  near  the  Oklahoma-Arkansas 
border  was  exactly  where  the  high  O3  on  the  first  flight  was  sampled, 
and  flight  time  was  very  close  to  the  time  of  the  satellite  observa¬ 
tion.  A  simple  calculation  using  the  ideal  gas  law  explains  that  a  40 
ppb  increase  in  Os  concentration  (which  was  measured  as  the  Sabreliner 
entered  the  dirty  air),  if  the  layer  were  4  to  6  km  thick,  would 
account  for  an  increase  of  15  Du  in  TCos.  In  addition,  the  small  spot 
of  low  TCoa  in  east  central  Oklahoma  may  be  evidence  of  the  gradient 
measured  later  in  the  day  at  10.6  Km.  The  picture  shows  vividly  the 
large  increase  in  TC03  north  of  the  jet  stream  resulting  from  the 
lower  tropopause  height. 

The  TOMS  sensor  measures  ultraviolet  (UV)  radiation  and  thus 
cannot  see  below  clouds.  Before  any  conclusion  could  be  drawn,  the 
effect  of  clouds  had  to  be  considered.  Comparison  between  the  TOMS 
imagery  and  infrared  (IR)  imagery  from  the  NOAA  GOES  satellite  for 
nearly  the  same  time  showed  not  only  that  the  maximum  in  TCo3  dis¬ 
cussed  above  was  near  the  clear  air  sampled,  but  that  wherever  there 
was  strong  convection,  TOMS  imagery  gave  lower  TC03.  NASA  uses  a 
simple  climatological  algorithm  to  correct  for  clouds  in  the  contin¬ 
uous  TOMS  database.  Was  the  variation  really  so  well  correlated  to 
the  tropospheric  features,  and  a  result  of  convection  pushing  up  the 
tropopause,  or  was  it  merely  an  anomaly  from  the  cloud  correction 
routine? 

The  algorithm  assumes  cloud  height  is  300  mb  in  the  tropics  and 
600  mb  at  the  poles;  at  intermediate  latitudes,  cloud  height  (in  mb) 
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is  given  as: 

300  +  (1  -  cos (2xlatitude)  x  150).  (1) 

The  lowest  700  mb  of  the  troposphere  is  assumed  to  contribute  19  Du 
(from  climatology) ,  distributed  linearly  with  pressure  (for  a  cloud  at 
300  mb  one  would  add  19  Du,  for  a  cloud  at  650  mb  one  would  add  9.5 
Du) .  If  ultraviolet  (UV)  reflectance  measured  by  the  sensor  is  less 
than  20  percent,  it  assumes  clear  skies;  greater  than  60  percent,  it 
assumes  a  cloud  at  the  climatological  height;  and  in  between  20  and  60 
percent,  the  correction  is  linearly  adjusted  between  zero  and  the 
total  adjustment  (P.K.  Bartia,  personal  communication,  1988) .  TOMS 
data  for  some  periods  is  available  from  NASA  corrected  using  measured 
IR  cloud  heights  and  climatological  tropospheric  O3 . 

Taking  the  grid  of  climatologically-corrected  data  (provided  by 
J.  Fishman,  personal  communication,  1988),  I  used  the  UV  reflectance 
from  each  pixel  to  subtract  the  NASA  correction  and  retrieve  the  raw 
TCos.  Horizontal  resolution  is  coarse,  with  just  8  pixels  in  every  10 
degrees  longitude,  and  1  pixel  for  every  degree  latitude.  The  raw 
data  was  then  corrected  with  a  new  algorithm. 

Our  algorithm  uses  better  estimates  of  tropopause  height,  cloud 
height  and  tropospheric  O3 .  The  Little  Rock  and  Oklahoma  City 
soundings  showed  the  tropopause  at  13  km.  Comparing  the  coldest 
temperatures  on  the  GOES  imagery  with  the  soundings  showed  the 
convection  close  to,  but  not  penetrating,  the  tropopause.  The 
strongest  UV  reflectance  closely  matched  the  coldest  IR  cloud  tops. 
For  the  strongest  convection  at  the  time  of  the  TOMS  data,  UV  reflec- 
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tance  was  88  percent.  Therefore,  an  altitude  of  13  km  was  set  to 
correspond  to  90  percent  reflectance,  and  0  km  (no  cloud)  was  set  to 
correspond  to  20  percent  reflectance.  Average  tropospheric  column  Oa 
was  calculated  by  averaging  the  three  best  fit  lines  from  the  O3 
profile  correlations;  two  lines  had  a  small  positive  slope  and  one  a 
small  negative  slope.  The  lines  were  extended  to  13  km  and  the  values 
at  each  0.1  km  were  averaged  over  the  three  profiles.  No  correction 
was  made  for  surface  elevation,  which  was  generally  below  1  km,  and 
the  expected  sharp  increase  in  O3  near  the  tropopause  was  not  measured 
or  included  (these  errors  should  partially  offset  each  other).  The 
resulting  average  was  49  ppb,  or  33.4  Du  from  the  surface  to  13  km. 
The  appropriate  amount  of  tropospheric  O3  was  added  to  the  raw  TCo3  by 
linear  interpolation  based  on  the  measured  reflectance.  No  attempt 
was  made  to  try  to  improve  on  a  linear  fit  as  reflectance  from  a  thin 
high  cloud  may  look  similar  to  a  solid  low  cloud,  with  the  former 
letting  through  some  radiation  from  below  and  the  latter  compl*  .y 
blocking  radiation  from  the  lowest  levels.  One  really  does  not  know 
what  levels  are  being  missed  by  the  TOMS  sensor. 

The  original  values  with  the  NASA  cloud  correction  and  the  newly 
corrected  TC03  values  are  shown  in  Figure  49.  The  new  algorithm 
increased  the  TC03  where  there  were  solid,  high  clouds.  While  high 
TC03  values  remain  on  the  Oklahoma-Arkansas  border,  it  is  no  longer  an 
isolated  maximum,  but  a  thin  band  of  elevated  TC03  extending  from  the 
high  TC03  to  the  northeast.  The  small  clean  spot  just  northwest  of 


this  band  remains. 


'yrv-'i’ 


INITIAL 


GRIDDED  TOMS  DATA 

(Unit  I)u  minus  200) 


89  90  91  &€  9 81  OBI  XUS3J3-1JB 

87  91  94  96  99i931^19?:UJLX27 

9  3  9  2 

9  2  9  3 
96  9  3  &$■& 

93  9i 

92  89  86  B 


CORRECTED 


89  90  91 
87  91  94 

93  92 


Hill 


5&€%0ifp 

mZMB 


■-99X&-3±'&-3 


Fig.  49.  Gridded  TOMS  data  from  22  June  1986.  The 
initial  data  was  derived  using  the  NASA  standard  correc¬ 
tion  for  tropospheric  O3 .  For  the  corrected  data,  the 
NASA  standard  correction  was  removed  and  a  new  correction 
was  made  based  on  the  observed  tropospheric  conditions  as 
described  in  the  text.  The  box  drawn  in  Figure  46 
indicates  the  geographical  extent  of  the  gridded  data. 
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While  the  corrected  TOMS  data  do  not  contradict  the  ideas 
proposed  about  the  origin  of  the  tropospheric  Oa  gradients,  they  fail 
to  provide  any  conclusive  support,  even  when  considered  along  with 
aircraft  data,  synoptic  data,  and  IR  (GOES)  satellite  imagery.  The 
TOMS  and  synoptic  data  are  on  a  very  coarse  scale  compared  to  the 
small  scale  features  sampled  by  the  Sabreliner.  The  clean  spot  in 
east  central  Oklahoma  is  only  one  pixel,  and  the  dirty  area  on  the 
Oklahoma-Arkansas  border  is  only  2  pixels;  the  synoptic  network 
consists  of  at  most  a  few  stations  per  state.  Upper  air  synoptic  data 
are  taken  twice  a  day;  TOMS  data  are  available  only  once  per  day,  and 
not  at  synoptic  observation  times.  In  addition,  the  TOMS  instrument 
measured  UV  reflectance  greater  than  20  percent  wherever  there  was  a 
hint  of  cloud  on  the  IR  imagery;  the  tropospheric  part  of  the  TC03  we 
are  comparing  with  the  aircraft  observations  was  often  partially 
blocked  by  clouds  and  corrections  based  on  average  tropospheric  Os 
further  obscure  real  tropospheric  gradients. 

The  coarse  time  resolution  of  the  TOMS  data  makes  interpretation 
of  the  movement  of  the  features  between  Figures  47  and  48  ambiguous. 
This  ambiguity  is  not  cleared  up  by  the  synoptic  data  close  to  the 
time  midway  between  the  two  TOMS  analyses.  In  the  equivalent  barotro- 
pic  assumption,  which  assumes  a  linear  increase  in  winds  from  the 
surface  to  the  tropopause,  synoptic  patterns  are  assumed  to  propagate 
at  the  mean  tropospheric  flow,  estimated  by  the  500  mb  winds.  But  in 
this  case  the  flow  is  strongly  non-barotropic  in  the  vicinity  of 
Oklahoma,  Arkansas  and  Texas  because  of  the  vertical  wind  shear 
(surface  winds  from  the  south  and  southeast  shifted  to  west  and  north- 
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west  by  300  mb).  In  this  case  waves  may  move  faster  than  the  parcels, 
or  remain  stationary  with  many  parcels  flowing  through.  Assuming  the 
feature  extending  from  the  Arkansas-Oklahoma  border  southwest  to  east 
central  Texas  in  Figure  48  is  associated  with  a  trough  in  the  tropo¬ 
sphere,  causing  a  depressed  tropopause  as  evidenced  by  the  high  TCoa, 
it  could  correspond  to  several  features  in  Figure  47,  24  hours 
earlier.  The  trough  could  be  the  feature  in  Figure  47  extending  from 
western  Arkansas  to  central  Texas,  which  remained  stationary  and 
deepened  slightly.  Or,  the  trough  could  have  pushed  south  from 
northern  Kansas.  Finally,  the  trough  could  have  broken  off  from  the 
larger  scale  trough  in  Figure  47,  which  at  that  time  extended  south  to 
northern  Colorado.  GOES  satellite  imagery  confirms  that  this  larger 
trough  pushed  south  to  central  Oklahoma,  and  is  seen  receding  north¬ 
ward  in  Figure  48.  The  wave  movement  occurs  independently  from  the 
individual  air  parcels. 


Analysis  of  all  the  available  data  suggests  the  following 
possible  interpretation  of  the  trace  gas  observations.  At  1200  GMT  on 
22  June,  the  upper  level  trough  was  moving  out  to  the  northeast. 
Winds  over  Oklahoma  were  light,  cutoff  from  the  faster  flow  to  the 
north  (Figure  1).  The  air  at  5.5  km  was  subsiding,  as  evidenced  by 
its  vertical  motion  and  large  dew  point  depression.  Although  the 
Sabreliner  flew  near  the  band  of  high  TCoa  through  central  Arkansas, 
it  only  sampled  high  local  Oa  in  the  small  area  near  the  Oklahoma 
border.  The  tropopause  was  depressed,  and  the  tropospheric  air  below 
was  sinking  along  this  band.  But  the  entire  band  was  not  pushed  down 
uniformly;  some  mesoscale  forcing  pushed  down  air  which  happened  to  be 
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outgassed  from  a  thunderstorm  (as  evidenced  by  the  high  CO  concentra¬ 
tion)  . 

As  for  the  large  gradient  sampled  at  10.6  km,  TOMS  data  confirms 
the  difference  in  TC03  at  the  origins  of  the  two  parcels  (Figures  46 
and  47).  The  air  at  10.6  km  moved  fast  enough  to  have  originated  from 
behind  the  convection  that  affected  the  region  earlier  in  the  day, 
while  the  air  slightly  lower  was  involved  in  the  recent  convection. 
The  fine  structure  at  10.6  km  was  caused  by  a  dynamic  process  which 
caused  confluence  of  the  parcels;  the  parcels  then  became  trapped  in 
the  slow  moving  flow  over  Oklahoma. 


V.  TROPOSPHERIC  OZONE  PRODUCTION 


Results  from  the  two  flights  on  22  June  suggest  some  of  the  O3 
sampled  is  of  tropospheric  origin,  either  produced  in  the  PBL  and 
brought  by  convection  into  the  free  troposphere,  or  produced  in-situ 
in  the  free  troposphere  from  ozone  precursors  (CO  or  hydrocarbons  in 
the  presence  of  NOx  and  sunlight).  This  section  estimates  the  tropo¬ 
spheric  O3  production  with  a  simple  photochemical  model. 

A  rough  estimation  can  be  made  of  whether  net  O3  production  is 
occurring  by  assuming  the  following  simple  reaction  scheme  (rate 
constants  from  NASA,  1985) : 

(4)  03  +  NO  -->  NOa  +  Oa  k< »1 . 8  (10- 1  * )  exp (-1370/T) 

(5)  NOa  +  hv  +  Oa  —  >  NO  +  O3  jo 

(7)  CO  +  OH  +  Oa  — >  COa  +  HOa  kT-1.5(10-‘3 ) (1+0.6P) 

(8)  HOa  +  NO  -->  NOa  +  OH  ka-3 .7 (10"1 2 ) exp(240/T) 

(9)  03  +  HOa  — >  OH  +  20a  ka»l .4 (lO"1 « ) exp(-580/T) 

The  change  in  O3  with  time  is  found  from  reactions  4,  5,  and  9: 


d[Oa]/dt  «  jo  [NOa]  -  k<[03][N0]  -  katOaJCHOa]  (2) 

And  the  change  in  NOa  with  time  is  found  from  reactions  4,  5,  and  8: 

d[NOa]/dt  -  ka  [NO]  [HOa ]  +  k«[NO][03]  -  jo  [NOa]  (3) 

Assuming  steady  state  in  NOa  yields: 

[NOalss  -  (ka [NO] [HOa]  +  ka [NO] [O3] ) / jo  (4) 

This  simplifies  equation  2  above: 

d[0a  ]  /dt  -  [HOa]  x  (ka  [NO]  -  kafOaj)  (5) 


If  the  expression  in  parentheses  is  positive,  O3  will  be  produced,  and 
if  negative  O3  will  be  destroyed.  This  expression  is  a  function  of 


temperature,  NO  concentration  and  Oa  concentration,  and  explains  bow 
one  can  define  a  critical  NO  concentration  for  O3  production  (e.g. 
Ridley  et  al.,  1987). 

To  see  if  Oa  production  was  likely  in  the  regions  sampled  on  the 
two  flights,  calculations  were  done  for  the  altitudes  of  level  flight 
where  NO  was  measured  (Table  3).  All  concentration  and  temperature 
data  at  a  level  were  averaged.  At  all  levels  where  NO  data  existed, 
Oa  production  was  possible.  In  most  cases,  as  altitude  increased,  Oa 
production  increased. 

Liu  et  al.  (1987)  estimated  the  total  column  production  of  O3 
over  the  United  States  in  summer  from  anthropogenic  NOx  to  be  about 
1(1012)  cm-2  s-1;  somewhat  greater  than  that  estimated  from  natural 
NOx,  which  they  estimate  to  range  from  1.9(10“)  to  12(10“)  cm-2  s_l. 
They  made  the  estimate  using  surface  measurements  from  the  rural  site 
at  Niwot  Ridge,  Colorado.  Average  NNHC  concentrations  were  derived 
from  measurements  by  Greenberg  and  Zimmerman  (1984).  Above  the  PBL, 
trace  gas  concentrations  were  set  at  values  representative  of  clean 
continental  air,  and  NNHC  concentrations  were  considered  negligible. 
To  compute  a  similar  estimation  with  this  data  set  required  a  more 
detailed  model  than  reactions  4-9,  to  consider  oxidation  of  both  CO 
and  hydrocarbons.  What  follows  is  an  order  of  magnitude  estimation  of 
the  column  Oa  production  under  the  conditions  found  on  each  flight, 
and  an  explanation  of  the  approximations  made. 

Table  4  lists  26  reactions  involved  in  tropospheric  O3  production 


TABLE  3 


OZONE  PRODUCTION  AT  ALTITUDES  OF  LEVEL  FLIGHT 

d03  »  [HO*]  (ks  [NO]  -  k»  [Os ] )  »  [HO*]  X 
dt 

If  X  >  0,  Os  produced.  If  X  <  0,  Os  destroyed 


Flight  1 

Altitude 

(m) 

X  (xlO-*) 

8800 

0.95 

10000 

4.8 

10600 

6.9 

9400 

2.4 

7900 

2.0 

5500 

0.44 

Flight  2 

Altitude 

(m) 

X  (xlO-*) 

7000 

0.63 

8800 

0.87 

9400 

0.98 

10600 

0.82 

ke  -  3.7(10-1*)  exp(240/T)  ;  ks  -  l.l(10-*«)  exp(-580/T) 


TABLE  4 


REACTIONS  INVOLVED  IN  OZONE  PRODUCTION 


(1)  HOt  +  NO  — >  NOs  +  OH 

(2)  Os  +  HOa  — >  OH  +  20* 

(3)  OH  +  0*  —  >  HO*  +  0* 

(4)  CHaO*  +  NO  -->  NO*  +  CHaO 

(5)  OOD)  +  HsO  —  >  20H 

(6)  CO  +  OH  +  Oa  -->  CO*  +  HO* 

(7)  0(1D)  +  0*  -->  0*  +  0(3P) 

(8)  OOD)  +  N*  —  >  Na  +  0(3P) 

(9)  Oa  +  hv  -->  0(*D)  +  0* 

(10)  CH4  +  OH  +  0*  -->  HaO  +  CHaO* 

(11)  CHaO  +  0*  -->  HaCO  +  HOa 

(12)  HaCO  +  OH  +  Oa  — >  HaO  +  HOa  +  CO 

(13)  HaCO  +  hv  -->  1/2 (2H0a  +  Ha )  +  CO 

(14)  HOa  +  HOa  — >  HaOa  +  Oa 

HOa  +  HOa  +  M  — >  HaOa  +  Oa  +  M 

(15)  HaOa  +  hv  -->  20H 

(16)  NOa  +  hv  — >  NO  +  0(3P) 

(17)  Oa  +  NO  — >  NOa  +  Oa 

(18)  0(3P)  +  Oa  -->  Oa 

(19)  OH  +  HOa  — >  HaO  +0a 

(20)  HOa  +  NOa  +  M  — >  HOaNO*  +  M 

(21)  HOaNO*  +  M  ~ >  HO*  +  NOa  +  N 

(22)  OH  +  NO*  +  M  — >  HNOa  +  M 

(23)  CHaO*  +  HO*  -->  CHaOOH  +  Oa 

(24)  CHaOOH  +  1/20H  -->  HaO 
+  1/2H*C0  +  l/2CHaO* 

(25)  HaOa  +  OH  -->  HaO  +  HO* 

(26)  HaOa  — >  rainout 


ki«3.7(10-»*)exp(240/T) 
k* -1 . 4 ( 10- »« ) exp (-580/T) 
ka  *1 . 6 ( 10“ 1  * ) exp ( -940/T ) 
k< *4.2(10“13)exp(180/T) 
ka“2.2(10~1#) 
k«-1.5(10“»3) (1+0.6P) 
kT»3.2(10“l 1 )exp(67/T) 
k* -1.8 (lOl-i »)exp(107/T) 

3* 

ki 0  *2 . 4 ( 10“ 1 3 ) exp (-1710/T) 
k»  i  -8 .4 (10- 1  * ) exp (-1200/T) 
kn*l. O(lO-n) 
jia-8.5(10-») 

ki<3>2.3  (10“ 13  )exp(590/T)  + 
1.7 (10" 3 3 )exp(1000/T)xM 
jia«l. 0(10-3) 

ki?*1.8(10-13)exp(-1370/T) 
ki a“6.0 (10“3 ) (T/300)-2 • 7xM 
kifl. 0(10-i‘)  (7+4p) 
k*o«2.0(10-3i ) (T/300)-2 • 7xM 
kai«kao/ (2.33 (10“27 )x 
exp(10870/T) ) 

k*a-2.6(10-30) (T/300)-3- 2xM 
kaa=7.7(10-i«)exp(1300/T) 
k*4»l. 0(10-i‘) 

ka b=3.1 (10-i 2 ) exp(1300/T) 
k* *=4. 0(10” 7 ) 


Source  of  rate  constants  is  NASA  (1985),  except  for  reaction  26,  which 
is  estimated  from  a  lifetime  of  one  month. 


and  destruction.  Pickering  (1987)  used  a  subset  of  the  most  important 
reactions  (reactions  1-10  on  Table  4) ,  to  derive  an  analytical 
solution  to  the  Os  production  problem.  The  full  set  of  reactions 
greatly  complicates  the  steady  state  analysis  that  follows,  and 
requires  a  numerical  solution.  Calculations  were  first  done  using 
the  analytical  model  below.  Then  some  results  were  compared  to  a 
complete  numerical  solution  using  a  model  recently  developed  by 
Pickering  and  Dickerson  (1988,  personal  communication).  (Note, 
reaction  numbers  refer  to  those  in  Table  4,  not  the  previous  text.) 
Os  concentration  changes  according  to  the  following: 

d[Os]/dt  *  [NO]  ([HOslki  +  [RO2  ]  k<  )  -  [00  D)  ]  [HaO]ko 

-  [Os]  ( [HOs ] ks  +  [OH] ks  )  (6) 


HOa  concentration  can  be  estimated  from: 


[HOi]  -  ktfeonOH]  +  2kiof  CH4]  [OH] 
ki  [NO]  +  ka[0s] 


(7) 


And  the  concentration  of  ROa  can  be  approximated  by: 


[CHsOa]  -  kio [CH« 1 [OH] 
k«  [NO] 


(8) 


Also,  OH  concentration  is  estimated  from  the  following: 


[OH]  -  2ka  [OOP)  1  fHaOl 
k«  [CO] 


(9) 


It  remains  to  estimate  the  concentration  of  00 D)  by: 


[OOD)] 


-i»IOs]_ 


(10) 


k7[0a]  +  ka[Na]  +  ks  [HaO] 

All  that  remains  unknown  is  ja ,  the  photolysis  rate  of  O3  by  radiation 
with  wavelength  less  than  318  nm.  This  was  measured  and  discussed  by 
Dickerson  et  al.  (1982).  This  quantity  is  highly  variable,  and  is  a 
function  primarily  of  altitude,  temperature,  pressure,  solar  angle 
(and  therefore  time  of  day  and  year),  aerosol  loading,  cloudiness, 
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total  column  ozone  and  surface  albedo.  Total  column  ozone  was  about 
300  Du  (from  TONS  data) .  Tropospheric  ozone  was  assumed  to  be  10 
percent,  and  the  effective  ozone  column  density  was  calculated  by 
dividing  300  Du  minus  a  fraction  of  10  percent,  based  on  altitude,  by 
the  cosine  of  the  solar  zenith  angle  for  the  aircraft  location  and 
flight  time.  Figure  4  in  Dickerson  et  al.  (1982),  a  graph  of  effec¬ 
tive  ozone  column  density  versus  photolysis  frequency  for  Boulder, 
Colorado  in  summer,  shows  that  measurements  and  calculations  from 
several  sources  fall  close  to  a  narrow  band.  This  figure  was  assumed 
to  be  valid  for  this  flight  location  as  well.  The  photolysis  fre¬ 
quency  from  this  curve  was  increased  to  include  the  effect  of  up- 
welling  radiation  (or  surface  albedo)  by  multiplying  by  the  net 
(upwelling  plus  downwelling)  radiation  and  dividing  by  downwelling 
radiation  (both  components  were  measured  on  the  Sabreliner  with  two 
Epply  UV  radiometers) .  The  inclusion  of  upwelling  radiation  increased 
j»  by  about  20  percent  in  most  cases,  but  went  as  high  as  a  55  percent 
increase  at  high  altitude  on  the  first  flight. 

Once  [0(1D)]  is  estimated,  one  can  estimate  [OH],  then  [HO2]  and 
[ROz ] ,  and  calculate  the  net  production  of  O3. 

The  atmosphere  was  divided  into  layers  to  perform  the  vertical 
integration  (see  Figures  10  and  21).  Layers  were  used  because 
continuous  data  did  not  exist  for  all  the  species  needed  in  the 
reaction  scheme  and  estimates  for  smaller  layers  did  not  seem  valid. 
This  introduces  error  into  the  model  as  the  product  of  two  means  is 
generally  not  the  same  as  the  mean  of  the  products  (Chameides  et  al., 


1987) .  All  data  measured  within  a  layer  were  averaged  without  regard 
to  exact  altitude  to  obtain  layer  averages  of  temperature,  pressure, 
CO,  03  and  H«0  concentrations.  For  example,  where  many  data  in  a 
layer  were  collected  at  a  constant  altitude,  the  variables  may  be 
biased  to  that  level  and  are  not  corrected  to  better  represent  the 
center  of  the  layer.  When  NO  was  not  measured  at  all  in  a  layer,  IS 
percent  of  NO*  concentration  was  used.  Liu  et  al.  (1987)  used  a 
(NO] / [NO* ]  ratio  of  1/3,  based  on  measurements  made  by  Williams  et  al. 
(1984)  at  Niwot  Ridge.  Their  instrument  did  not  have  a  PAN  interfer¬ 
ence,  nonetheless  15  percent  is  a  conservative  estimate.  This 
estimate  inputs  substantial  error  into  the  calculation  of  03  produc¬ 
tion,  as  the  production  of  O3  goes  up  with  the  ratio  NO2  to  NO.  The 
tropospheric  column  production  is  an  underestimate,  since  the  tropo¬ 
sphere  above  10.5  Am  was  not  included  in  the  calculation  as  it  was  not 
sampled.  This  is  partially  balanced  by  the  PBL  values,  which  may  be 
too  high,  as  no  attempt  was  made  to  include  surface  loss  of  O3 . 

Although  NMHC  concentrations  are  .sually  several  orders  of 
magnitude  lower  than  that  of  CH< ,  their  effect  on  atmospheric  chemis¬ 
try  is  important  because  they  react  much  more  quickly  than  CH4  (Green¬ 
berg  and  Zimmerman,  1984).  Unlike  the  remote  Pacific  (Chameides  et 
al.,  1987),  NMHC  concentrations  were  significant  in  the  free  tropo¬ 
sphere  of  the  central  US,  and  neglecting  NMHCs  can  underestimate  O3 
production  (Liu  et  al.,  1987).  Estimates  by  Zimmermann  et  al.  (1978) 
and  Hanst  et  al.  (1980),  indicate  the  atmospheric  breakdown  of  NMHCs 
from  vegetation  is  a  major  source  of  atmospheric  CO.  No  grab  samples 
of  hydrocarbons  were  collected  during  either  flight  on  22  June  1985, 
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but  throughout  the  PRE-STORM  Project  twelve  different  samples  were 
collected  under  various  conditions  and  analyzed  by  J.P.  Greenberg  and 
P.R.  Zimmerman  of  NCAR.  The  samples  were  analyzed  using  gas  chroma¬ 
tography  with  a  flame  ionization  detector,  as  discussed  in  Greenberg 
and  Zimmerman  (1984) .  Basically  the  technique  consists  of  passing  the 
sample  through  various  analytical  columns  along  with  a  carrier  gas. 
Each  hydrocarbon  moves  through  the  column  at  a  different  rate,  and 
will  show  up  as  a  peak  on  a  chromatogram  at  a  representative  retention 
time.  The  mixing  ratio  can  be  inferred  from  the  height  of  the  peak  or 
the  area  under  the  peak. 


Using  the  measured  hydrocarbon  concentrations,  the  following 
simple  scheme  was  used  to  approximate  the  effective  methane  concentra¬ 
tion: 

[CIU ]• f i  -  [CH«]  +  £  (n[CnH]  X  kca/kci<)  (11) 

Where  n[CaH]  indicates  the  concentration  of  a  specific  hydrocarbon 
multiplied  by  the  number  of  carbon  atoms  it  contains,  ken  is  the  rate 
of  reaction  of  that  hydrocarbon  with  OH,  and  kca<  is  the  rate  of 
reaction  of  CH«  with  OH.  Table  5  lists  the  number  of  times  each 
hydrocarbon  was  found  in  a  sample  and  the  rate  constants  used.  The 
gas  chromatography  technique  fails  to  detect  NMHCs  with  more  than  10 
carbon  atoms  or  oxygenated  hydrocarbons  (Liu  et  al.,  1987)  The 
assumption  is  that  the  hydrogen  atom  extraction  is  the  rate  determin¬ 
ing  step,  and  after  that  the  rest  of  the  molecule  reacts  almost 
instantaneously.  Analysis  of  the  initial  results  found  that  where 
propylene  was  detected,  it  had  an  extremely  large  effect  on  this 
calculation.  The  propylene  reaction  with  the  first  OH  is  extremely 
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TABLE  5 


HYDROCARBON  ANALYSIS 


SPECIES  NUMBER  OF  SAMPLES 

RATE  CONSTANT 

Methane 

12 

8.0(10-19)exp(-1710/T)  k 

CO 

12 

n/a 

Ethane 

12 

1.1 (10~ 11 )exp(-1090/T)  n 

Ethylene 

9 

8 .4 (10_1 2 ) exp (8 . 39/T)  k 

Acetylene 

10 

2.0 (10_ 12 ) exp(-250/T)  h 

Propane 

12 

1.6 ( 10" 2 1 ) exp(-800/T)  n 

Propylene 

8 

see  below  e,k 

Isobutane 

11 

2.4(10-12)exp(-390/T)  k 

N-Butane 

12 

2.6 (10- 12 )exp(-560/T)  k 

Isopentane 

10 

4 . 1 (10" 1 2 ) exp (-530/T)  e 

Neopentane 

11 

4.1 ( 10" 1 2 ) exp (-530/T)  e 

N-Pentane 

12 

4.1 (10"12 )exp(-530/T)  k 

Isoprene 

2 

9.9(10-n)exp(500/T)  k 

1-Pentene 

2 

10" 10  e 

Neohexane 

10 

5.7 (10“12 )exp(-530/T)  e 

4-Me-l-Pentene 

4 

10- 10  e 

2-Me-Pentane 

10 

5.7 (10“12 ) exp(-530/T)  e 

3-Me-Pentane 

9 

5.7 (10“12 ) exp (-530/T)  e 

1-Hexene 

3 

10"10  e 

N-Hexane 

11 

5.7(10"12 )exp(-530/T)  k 

Benzene 

12 

7.57 (10“12 )exp(-529/T)  p 

Cyclohexane 

9 

5.7 (10"12 ) exp(-530/T)  e 

2,4Di-Me-Pentane 

1 

7.3 (10“12 )exp(-530/T)  e 

3-Me-Hexane 

1 

7.3 (10" 1 2 ) exp(-530/T)  e 

2 , 3-Di-Me-Pentane 

2 

7.3 ( 10" 12 ) exp(-530/T)  e 

1-Heptcne 

3 

10"10  e 

N-Heptane 

9 

7.3 (10"ia ) exp(-530/T)  k 

Me-Cyclohexane 

7 

7.3 (10-1 2 )exp(-530/T)  e 

Toluene 

10 

10-io  e 

2,2, 4Tri-Me-Pentane 

2 

9. 0(10" 12 )exp(-530/T)  e 

2,3, 4Tri-Me-Pentane 

1 

9.0 (10- 1 2 ) exp (-530/T)  e 

3-Me-Heptane 

2 

9.0 (10"i 2 ) exp(-530/T)  e 

N-Octane 

8 

9.0(10-io)exp(-530/T)  k 

1-Octene 

6 

10-1°  e 

Ethyl-Benzene 

2 

10-n  e 

P-Xylene 

4 

2.4(10-ii)  k 

M-Xylene 

1 

2.4(10-ii)  k 

0-Xylene 

2.4(10-ii)  k 

1-Nonene 

7 

10-1°  e 

N-Nonane 

9 

10-n  e 

Isopropylbenzene 

2 

10-io  e 

N-Propyl-Benzene 

3 

10-1°  e 

1,3,5-Tri-Me-Tol 

2 

10-n  e 

1-Decene 

1 

10-io  e 

N-Decane 

8 

10-n  e 

Alpha-Pinene 

1 

5.6(10-ii)  k 

Beta-Pinene 

2 

10-io  e 

Delta-3-Carene 

3 

10-io  e 

Delta-4-Carene  1 

Betaphellandrene  3 

Myrcene  1 

Butanal  6 

Pentanal  9 

Hexanal  6 

3-Heptanone  3 

Heptanal  8 

Octanal  2 

Methychloroform  12  ! 

e  -  estimated 

k  -  Kerr  and  Calvert  (1984) 
h  -  Hampson  (1980) 

p  -  Finlayson-Pitts  and  Pitts  (1986) 
n  -  NASA  (1985) 

(Rate  constant  for  propylene  used  :  : 
+  10-“  +  2.9(10-“  )exp(540/T)]) 


10-“  e 
10-“  e 
10-“  e 

2.5 (10" 1 1 ) exp (250/T)  k 

1.6(10-“  )exp(260/T)  k 

10-“  e 

10-“  e 

10-“  e 

10-“  e 

5.4(10-“)exp(-1820/T) 


1/3  [1.6 (10" 1 1 ) exp(260/T) 
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fast  (better  than  2(10~i0)  s- 1  at  rooa  tesperature  and  pressure).  It 
was  decided  to  slow  the  breakdown  rate  by  averaging  this  rate  with  the 
breakdown  rates  of  two  possible  stable  intermediates;  this  is  the  rate 
given  for  propylene  as  a  footnote  to  Table  5. 

The  total  NMHC  concentrations  and  the  ratios  of  effective  methane 
to  measured  methane  concentration  are  listed  for  each  sample  in  Table 
6.  The  smallest  ratio  was  1.8  for  a  very  clean  sample  and  the  largest 
was  a  tremendous  549  from  an  extremely  dirty  sample. 

Certainly  this  scheme  is  an  over  simplification.  Exact  calcula¬ 
tions  of  the  effect  of  NMHCs  would  require  a  detailed  scheme  of  the 
complete  breakdown  rate  of  each  hydrocarbon,  checking  especially  for 
stable  intermediates.  Also,  some  NMHCs  may  breakdown  to  CO2  without 
becoming  CO.  Nevertheless,  this  exercise  suggests  their  effect  is  far 
from  negligible  in  increasing  [CH4]«ff.  The  ratios  showed  tremendous 
variation,  but  were  not  correlated  with  altitude  or  CO  concentration. 
Despite  their  relatively  short  lifetimes,  NMBCs  were  observed  at  high 
enough  concentrations  in  the  free  troposphere  to  greatly  increase 
[CH4 ]• f f  in  most  samples.  Host  hydrocarbons  have  longer  lifetimes  at 
the  lower  temperatures  and  pressures  at  high  altitudes.  Some  past 
studies  only  included  the  effect  of  NMBCs  on  PBL  Os  production 
(Vukavich  et  al.,  1985;  Trainer  et  al.,  1987).  Data  suggest  that 
NMHCs  may  effect  03  production  throughout  the  free  troposphere  when 
they  escape  the  PBL  by  means  such  as  convection. 

Methane  concentrations,  because  of  its  long  atmospheric  lifetime, 


are  almost  always  close  to  1.6  ppm  in  the  troposphere.  For  the  Os 
column  production,  it  was  decided  to  try  three  effective  methane 
concentrations,  1.6  ppm  (ignoring  the  effect  of  NMBCs) ,  2.9  ppm  (a 
ratio  of  1.8,  equal  to  the  cleanest  case  measured)  and  120  ppm  (a 
ratio  of  75,  but  cleaner  than  the  three  dirtiest  samples  measured). 
The  same  was  used  in  all  four  layers,  and  three  calculations 
were  done  for  each  flight.  It  is  reasonable  to  expect  at  least 
moderate  values  of  man-made  NMHCs  in  all  layers  of  the  first  flight 
and  the  lowest  two  layers  of  the  second  flight;  low  level  wind  trajec¬ 
tories  were  from  the  south  through  southwest,  bringing  the  air  from 
the  Gulf  of  Mexico  through  east  Texas  and  Louisiana,  which  has  many 
refineries  and  power  plants.  As  reported  in  the  1980  NAPAP  (National 
Air  Pollution  Assessment  Program)  emission  inventories  by  Wagner  et 
al.  (1986),  this  region  has  high  emissions  of  NO*  and  SOz;  which 
suggests  accompanying  high  anthropogenic  hydrocarbon  emissions.  As 
discussed,  this  low  level  air  was  mixed  to  the  higher  layers  by 
Convection. 

In  all  six  calculations,  local  photochemistry  produced  O3  in 
every  layer.  Table  7  shows  the  net  production  rate  in  each  layer  and 
the  sum  for  the  column  in  each  case. 

Production  during  the  later  flight,  as  solar  zenith  angle  was 
rapidly  increasing,  is  less  than  1/7  of  that  during  the  first  flight, 
when  the  sun  was  nearly  overhead.  To  make  a  rough  estimate  of  the 
diurnal  average  production  rates,  the  calculations  were  redone  with 
zenith  angle  of  45  degrees,  and  the  results  were  divided  by  two.  The 


TABLE  7 


NET  OZONE  COLUMN  PRODUCTION 
Units  :  10ia  nolec/ (cm2s) 


Flight  1 

[CH« ] • f  r  (ppm) 

1.6 

2.9 

Layer  1: 

3.3 

5.4 

Layer  2: 

1.3 

2.2 

Layer  3: 

.75 

1.2 

Layer  4: 

.17 

.25 

Total: 

5.5 

9.1 

Diurnal  Average: 

1.8 

2.9 

Numerical  Model: 

1.0 

1.3 

Flight  2 

[CH« ] • t  f  (ppm) 

1.6 

2.9 

Layer  1: 

.47 

.76 

Layer  2: 

.19 

.31 

Layer  3: 

.052 

.086 

Layer  4: 

.024 

.037 

Total: 

.73 

1.2 

Diurnal  Average: 

1.4 

2.2 

Numerical  Model: 

1.3 

1.7 

correction  for  upwelling  radiation  was  not  changed,  which  assumes  the 
surface  albedo  is  independent  of  zenith  angle.  The  results  of  these 
calculations  are  also  presented  in  Table  7. 


To  check  the  accuracy  of  this  simple  model,  the  diurnal  averages 
were  also  calculated  using  the  numerical  model  which  includes  the  26 
reactions  in  Table  4.  The  model  numerically  solves  the  differential 
equations  necessary  to  arrive  at  the  O3  production.  One  must  make 
initial  guesses  of  OH,  H0»,  CHaO,  HiCO,  CHaOs,  0(3P),  HO2NO2,  and 
CHaOOH  concentration,  but  the  model  recalculates  these  quantities  and 
seems  to  be  stable  to  small  changes  of  these  input  parameters.  The 
greatest  source  of  uncertainty  in  the  model  comes  from  the  necessity 
to  guess  NOa  concentration  from  the  measurement  of  NO*  (this  is 
complicated  by  the  PAN  interference).  At  higher  altitudes,  more  odd 
nitrogen  would  be  expected  to  be  in  one  of  the  reservoir  species  such 
as  PAN  than  as  NOa ,  so  I  used  a  decreasing  percentage  with  altitude  of 
the  layer  average  NOx  to  estimate  NOa  (60  percent  in  layer  I,  50 
percent  in  layer  II,  40  percent  in  layer  III,  and  20  percent  in  layer 
IV).  The  photolysis  rate  of  NOa,  jie,  was  estimated  from  an  equation 
based  on  zenith  angle  in  Dickerson  et  al.  (1982). 


The  results  of  the  numerical  model  are  also  presented  in  Table  7. 
For  the  case  with  only  CH< ,  the  results  of  the  analytical  and  numeri¬ 
cal  solution  agree  quite  well.  At  the  higher  [CH4]«ff  ratios,  the 
numerical  solution  tends  to  moderate  the  effect  of  high  NMHCs.  It 
appears  high  concentrations  of  NMHCs  may  increase  net  column  03 
production  by  a  factor  of  4  or  5  over  that  expected  with  only  CH4 , 
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rather  than  by  more  than  an  order  of  magnitude  as  suggested  by  the 
analytical  model. 

The  production  estimates  using  the  lowest  two  [ClUlaff  ratios  are 
of  about  the  same  magnitude  as  the  estimate  given  above  by  Liu  et  al. 
(1987)  (sum  of  anthropogenic  and  natural  production) ,  who  assumed  no 
NMHCs  in  the  free  troposphere;  inclusion  of  high  NMHCs  increases  the 
column  production  by  a  factor  of  two  to  four  over  that  estimated  by 
Liu  et  al.  (1987).  These  results  are  two  orders  of  magnitude  larger 
than  the  Os  production  calculated  for  the  remote  Pacific  by  Chameides 
et  al.  (1987).  To  get  a  feel  for  the  importance  of  tropospheric 
production,  2 ( 101 2 )  molec/cmas  in  the  mid-troposphere  (600mb,  270K) 
corresponds  to  production  of  about  10  ppb/day.  The  average  tropo¬ 
spheric  Os  concentration  is  around  50  ppb,  so  for  the  day  sampled, 
tropospheric  production  of  Os  was  significant. 


VI.  SUMMARY  AND  CONCLUSIONS 


This  study  presents  some  of  the  first  aircraft  data  analyzed  from 
flights  over  non-urban  North  America  where  Os ,  CO  and  odd  nitrogen 
were  measured  concurrently  in  the  free  troposphere.  The  two  flights 
presented  complement  each  other;  the  first  shows  a  troposphere  well 
mixed  by  recent  convection,  and  the  second  shows  a  stable  stratifica¬ 
tion  several  hours  later.  The  flight  area  did  not  experience  a 
frontal  passage;  low-level  winds  were  from  the  south  through  south¬ 
west,  shifting  to  westerly  aloft  for  both  flights. 

Analysis  of  the  first  flight  showed  that  immediately  after 
intense  convection,  large  areas  (the  flight  traversed  more  than  500 
km)  of  the  atmosphere  are  disturbed  by  mesoscale  mixing.  The  tropo¬ 
sphere  is  not  thoroughly  mixed,  but  large  volumes  of  air  have  been 
displaced  vertically,  resulting  in  profiles  very  different  than  those 
predicted  by  eddy  diffusion.  When  there  is  a  large  amount  of  convec¬ 
tion,  horizontal  variations  will  be  large,  otherwise  vertical  varia¬ 
tions  will  predominate  as  on  the  second  flight.  The  effects  of 
thunderstorms  and  convection  should  not  be  ignored  in  models  of 
tropospheric  trace  gas  concentrations. 

The  depth  of  tropospheric  air  with  concentrations  typical  of  the 
PBL  can  be  very  deep  (greater  than  5  km) ,  during  summer  in  the  Central 
US.  Long-range  transport  is  increased  at  these  higher  altitudes,  and 
the  lifetimes  of  O3  precursors,  such  as  NOy ,  are  greater.  Like 
thunderstorms,  the  mechanism  that  causes  this  raised  PBL  has  potential 


to  increase  net  tropospheric  O3  production. 

A  strong  gradient  in  O3 ,  caused  by  neither  tropopause  fold  or  a 
thunderstorm,  was  observed  at  10.6  km  on  the  second  flight  in  eastern 
Oklahoma.  The  mesoscale  gradient  was  not  correlated  with  CO  or 
potential  temperature,  but  was  strongly  correlated  with  winds,  showing 
that  synoptic  scale  dynamic  phenomena  causing  confluence  of  air 
parcels  can  dramatically  alter  trace  gas  concentrations  in  the  upper 
free  troposphere. 

Tropical  tropospheric  O3  features  studied  by  Fishman  et  al. 
(1986)  were  often  the  dominant  features  in  the  TCo3.  But  in  this 
case,  analysis  of  TOMS  data,  after  correction  for  the  cloud  compensa¬ 
tion  algorithm,  failed  to  provide  evidence  to  either  support  or  refute 
whether  the  tropospheric  03  features  sampled  on  these  missions  caused 
the  gradients  in  TC03  seen  on  TOMS  imagery.  The  TOMS  imagery  can  be 
used  with  other  data  to  understand  the  large  scale  synoptic  features; 
however,  it  must  be  used  with  caution.  Over  continental  areas  in 
summer,  the  cloud  algorithm  tends  to  underestimate  the  contribution  of 
tropospheric  O3  by  neglecting  the  very  high  cloud  tops  from  thunder¬ 
storm  activity  and  underestimating  the  background  tropospheric  O3 
concentration.  As  a  result,  clear  areas  will  tend  to  show  up  as 
maxima,  and  cloudy  areas  as  minima,  where  no  real  gradient  exists. 
This  effect  should  be  considered  when  TOMS  data  are  used  for  applica¬ 
tions  other  than  estimating  total  global  column  Os  or  comparing  the 
dominant  stratospheric  features. 


A  simple  photochemical  model  found  in-situ  production  possible  at 
all  altitudes  under  the  conditions  measured  on  each  flight.  Net 
tropospheric  column  production  without  NMHCs  was  close  to  the  estimate 
of  Liu  et  al.  (1987).  Consideration  of  NHHCs  in  a  simplified  "effec¬ 
tive  CH«"  scheme  shows  that  NMHCs  have  an  effect  on  tropospheric  Os 
production  due  to  their  strong  reactivity,  and  increased  production  by 
a  factor  of  two  to  four  over  Liu  et  al.  (1987).  Tropospheric  produc¬ 
tion  appeared  to  be  an  important  source  of  tropospheric  O3  on  this 
day,  and  was  enhanced  by  the  strong  convective  mixing  of  Os  precursors 
and  very  deep  PBL.  NMHC  concentrations  were  not  correlated  with 
altitude  and  high  NMHC  concentrations  were  measured  above  the  PBL. 


APPENDIX  A.  DATA  PROCESSING 
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Computer  programs  described  by  Pickering  (1987)  were  used 
initially  to  process  the  minimally-processed  data  received  on  tape 
from  NCAR.  The  output  of  these  programs  was  reviewed  and  the  data 
quality  improved  as  described  below. 

1.  Nitrogen  Oxides  | 

The  raw  data  were  processed  by  subtracting  the  recorded  voltages 
from  the  baseline  made  by  linearly  interpolating  between  the  average 

( 

values  of  two  consecutive  backgrounds.  Outgassing  from  the  converters  | 

| 

was  checked  by  cycling  zero  air  through  the  channels.  Ideally,  any  j 

deviations  of  the  zero  air  from  the  background  should  be  subtracted 

from  the  appropriate  channel.  It  was  determined  that  the  zero  air 

values  were  negligible,  but  there  was  enough  noise  in  the  zero  air 

reading  to  cause  high  uncertainty  in  measurements  near  the  detection 

limit.  Due  to  the  linear  interpolation  and  the  background  noise,  some 

of  the  NO  values  were  calculated  to  be  negative  after  the  first 

processing.  When  this  happened,  the  data  were  analyzed  by  hand,  and 

corrections  were  made  to  account  for  any  erroneously  high  voltages  in 

the  background  or  what  appeared  to  be  an  obvious  non-linear  drop  in 

background  values  when  the  instrument  was  still  warming  up.  Positive, 

but  small  values  were  extracted  for  NO.  Several  strong  spikes  were 

eliminated  from  the  analyzed  data;  most  were  NOy ,  immediately  after  a 

mode  change.  One  very  sharp  spike  of  NO  (suspected  to  be  caused  by  a 

faulty  solenoid)  on  the  second  flight  at  1809  CDT  in  clear  air  was 

eliminated  as  well.  At  the  end  of  the  second  flight,  during  the 
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second  penetration  into  the  boundary  layer,  it  was  determined  some 
zero  air  was  getting  into  the  sample  flow.  The  measurements  were  not 
used  for  correlations  or  calculations.  Before  the  first  background 
and  after  the  last  background,  the  data  were  not  processed. 


2 .  Ozone 

Little  additional  processing  was  required  on  the  continuous  O3 
data  set.  Some  measurements  at  the  beginning  of  each  flight  were 
ignored  while  the  detector  was  warming  up.  Due  to  spiking  in  the  data 
acquisition  system,  several  points  on  each  flight  were  cut:  on  the 
first  flight  near  1717  and  1806  CDT,  and  on  the  second  flight  near 
1304  and  1454  CDT. 


3.  Carbon  Monoxide 

As  with  the  nitrogen  oxide  detector,  data  had  to  be  processed  by 
subtracting  the  voltages  from  the  baseline  drawn  between  the  average 
values  of  two  consecutive  backgrounds.  A  30  second  response  time  was 
assumed  when  processing  the  data  to  account  for  the  time  required  to 
pump  air  through  the  detector  (O3  and  odd  nitrogen  detection  are 
nearly  instantaneous).  This  lag  was  checked  against  O3  and  water 
vapor  in  cases  where  the  Sabreliner  was  entering  and  leaving  the 
boundary  layer  (ie.,  where  a  rapid  change  would  be  expected  in  each 
variable) .  The  changes  in  the  three  variables  showed  a  reasonable 
match. 
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For  the  first  flight,  data  were  eliminated  near  1828,  between 
1834  and  1842,  and  after  1851  CDT.  For  the  second  flight,  a  spike  in 
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one  background,  caused  by  a  sudden  change  in  aircraft  attitude 
(Dickerson  and  Delany,  1988),  was  replaced  with  a  more  reasonable 
value  so  CO  concentrations  would  be  positive. 


APPENDIX  B 


Average  data  for  every  1  minute.  Missing  data  is  indicated  by  0.0. 
The  NO  column  contains  NO  (status  1.),  NOx  (status  2.)  and  NOy 
(status  3 . )  . 

FLIGHT  1 


HR 

MIN 

ALT 

LAT 

LONG 

TEMP 

TD 

WD 

WS 

03 

CO 

NO  STAT 

(M)  . 

(DEG  N) 

(DEG  E) 

(C) 

(C) 

(DEG) 

(M/S) 

(PPB) 

(PPB) 

(PPT) 

12 

19 

363.1 

0.00 

0.00 

30 . 2 

17 . 7 

0  . 

0.0 

0.0 

0 . 0 

0  . 

0  . 

12 

20 

418.8 

0.00 

0 . 00 

27.8 

20.0 

210  . 

4.2 

0 . 0 

0 . 0 

0  . 

0  . 

12 

21 

990 . 5 

35.34 

-97.60 

21.8 

18 . 0 

189  . 

5 . 2 

21.5 

0 . 0 

0  . 

0  . 

12 

22 

1682 . 3 

35.33 

-97 . 54 

18 , 4 

10 .6 

169. 

5.8 

26 . 8 

0 . 0 

0  . 

0  . 

12 

23 

2263.9 

35.35 

-97.46 

15  9 

6 . 3 

142  . 

3.0 

32.6 

0.0 

0  . 

0  . 

12 

24 

2801.2 

35.38 

-97 . 37 

12.5 

1.7 

119. 

5 . 2 

33. 3 

0.0 

0  . 

0  . 

12 

25 

3485.3 

35.39 

-97 . 28 

8 . 2 

-5.7 

51  . 

4.5 

46 . 1 

0 . 0 

0  . 

0  . 

12 

26 

4057 . 5 

35.38 

-97 . 19 

2 . 8 

-7.9 

183  . 

5.6 

49 . 5 

0 . 0 

0  . 

0  . 

12 

27 

4568.4 

35.36 

-97 . 08 

-1.5 

-7.6 

330  . 

5.6 

46.4 

0 . 0 

0  . 

0  . 

12 

28 

4991.3 

35.34 

-96.97 

-3.3 

-14.2 

308  . 

5.7 

41 .4 

124.4 

0  . 

0  . 

12 

29 

5480 . 0 

35.34 

-96.88 

-5.0 

-22.0 

245. 

4 . 1 

31 . 1 

113.1 

0  . 

0  . 

12 

30 

5892.8 

35.34 

-96 . 76 

-7 . 7 

-25.1 

253  . 

4.8 

38 . 5 

98.8 

0  . 

0  . 

12 

31 

6140 . 3 

35.32 

-96.65 

-10.1 

-26.9 

325. 

5.6 

42 . 7 

93 . 5 

0  . 

0  . 

12 

32 

6480.0 

35.31 

-96.54 

-13.1 

-28.4 

293  . 

7 . 3 

52 . 0 

89.6 

0  . 

0  . 

12 

33 

6872.3 

35.30 

-96.44 

-16.3 

-26.2 

277  . 

8 . 2 

63.9 

91 . 2 

0  . 

0  . 

12 

34 

7184.5 

35.29 

-96.32 

-18.5 

-24.5 

297  . 

6.6 

46 . 2 

98 . 3 

0  . 

0  . 

12 

35 

7406.7 

35.28 

-96 . 19 

-20.1 

-25.4 

291  . 

7 . 0 

42.9 

99 . 2 

0  . 

0  . 

12 

36 

7603.2 

35.27 

-96.08 

-22.0 

-26.2 

284  . 

8 . 5 

47.6 

0 . 0 

0  . 

0  . 

12 

37 

7595.8 

35.23 

-95.96 

-22.1 

-26.0 

288  . 

8 . 3 

46 . 3 

0 . 0 

0  . 

0  . 

12 

38 

7590.5 

35.18 

-95..  82 

-21 . 8 

-26.2 

286. 

7. 1 

48.5 

0 . 0 

0  . 

0  . 

12 

39 

7840.7 

35.14 

-95.71 

-23.5 

-27.3 

294  . 

8.2 

48.8 

92.4 

0  . 

0  . 

12 

40 

8250.3 

35.11 

-95.60 

-26.5 

-29.0 

279  . 

7 . 0 

51 . 0 

124 . 2 

0  . 

0  . 

12 

41 

8502 . 3 

35.10 

-95.46 

-27.4 

-31.9 

261  . 

8.5 

56 . 0 

121 . 3 

0  . 

0  . 

12 

42 

8663.4 

35 . 09 

-95.32 

-27.8 

-35.2 

254  . 

11.9 

53.2 

120 . 2 

0  . 

0  . 

12 

43 

8829.8 

35.08 

-95.21 

-29.1 

-38.0 

254. 

12.8 

48.9 

125.4 

0  . 

0  . 

12 

44 

8804.6 

35.06 

-95.06 

-29.0 

-39.3 

253. 

12.4 

52 . 3 

128.4 

0  . 

0  . 

12 

45 

8810.6 

35.04 

-94.92 

-29.3 

-37.3 

259. 

11 . 1 

54.8 

116 . 3 

0. 

0  . 

12 

46 

8802.0 

35.03 

-94.79 

-29.4 

-36.3 

267. 

9.9 

51.5 

124.7 

1964. 

3. 

12 

47 

8796.7 

35.00 

-94.65 

-29.2 

-36.3 

274  . 

8.7 

48.5 

121 . 1 

1497. 

3. 

12 

48 

8787 . 3 

34.98 

-94.50 

-29.1 

-35.4 

276  . 

7.6 

45.6 

0.0 

1612  . 

3  . 

12 

49 

8787.4 

34.97 

-94.35 

-29.3 

-36.3 

272  . 

6.3 

43.8 

0 .0 

0. 

0  . 

12 

50 

8794.0 

34.95 

-94.22 

-29.4 

-36.5 

275  . 

7.6 

47.0 

0 . 0 

0  . 

0  . 

12 

51 

8791.9 

34.93 

-94.08 

-29.4 

-35.0 

282  . 

6.7 

45.0 

115.5 

0  . 

0  . 

12 

52 

8792.9 

34.91 

-93.93 

-29.4 

-35.6 

290  . 

5.4 

45.3 

92.6 

87. 

1  . 

12 

53 

8792.3 

34.90 

-93.81 

-29.5 

-34.5 

283. 

5.4 

45.8 

96 . 2 

9i  . 

1  . 

12 

54 

8791.4 

34.88 

-93.68 

-29 . 5 

-34.2 

253. 

6.6 

41 . 5 

93.6 

155  . 

1  . 

12 

55 

8789.9 

34.86 

-93.53 

-29.7 

-33.7 

241  . 

6 . 1 

37.9 

94.2 

517  . 

2  . 

12 

56 

8796.6 

34.85 

-93.40 

-29.9 

-33.9 

244  . 

8.3 

38.4 

92.4 

626  . 

2  . 

12 

57 

8804.1 

34.83 

-93.26 

-30.0 

-33.5 

243  . 

9.7 

39.6 

88 . 0 

1064  . 

2  . 

12 

58 

9047.4 

34.82 

-93.13 

-30.9 

-34.3 

265. 

3  0.4 

39.6 

96.6 

1593. 

3. 

12 

59 

9460.7 

34.81 

-92.99 

-33.9 

-39.1 

269. 

15.8 

38.4 

104.4 

1281  . 

3  . 

13 

0 

9752.6 

34.79 

-92.88 

-36.6 

-42 . 5 

272  . 

19.9 

43 . 1 

0 . 0 

1174. 

3. 

13 

1 

9955.0 

34.78 

-92.78 

-38.9 

-42.9 

266  . 

20.6 

39.5 

0.0 

0  . 

0  . 

13 

2 

10040.6 

34.77 

-92.65 

-39.9 

-43.8 

268. 

21.6 

39.2 

0 . 0 

0  . 

0  . 

13 

3 

10024.4 

34.75 

-92 . 53 

-39.8 

-44.5 

270  . 

15 . 2 

3  1  2 

102.5 

0  . 

0  . 

13 

4 

10029.0 

34.74 

-92.42 

-39.8 

-44 . 1 

265  . 

17 . 3 

38.3 

111.6 

588. 

1  . 

13 

5 

10021.5 

34.72 

-92 . 28 

-39.9 

-43 . 2 

262  . 

17 . 8 

35.7 

104.6 

688. 

1  . 

13 

6 

10012.6 

34.70 

-92.14 

-40 . 3 

-42.8 

271  . 

12.8 

39 . 2 

101 . 3 

305. 

1  . 

13 

7 

10004.6 

34.67 

-92.01 

-40.4 

-43 . 3 

282  . 

10.4 

49.3 

91.7 

548. 

2  . 

13 

8 

10274.0 

34.64 

-91.89 

-41.8 

-45.7 

302  . 

17.7 

48.5 

115.8 

480  . 

2  . 

13 

9 

10611.4 

34.61 

-91.76 

-44.6 

-49 . 8 

293. 

24.7 

44 . 3 

116.2 

651  . 

2  . 

13 

10 

10634.4 

34.59 

-91.63 

-45.1 

-50 . 9 

279. 

23.4 

46.5 

102.6 

1754  . 

3  . 

13 

11 

10632.9 

34.58 

-91 . 50 

-45.2 

-52 . 2 

277  . 

18 . 5 

42 . 8 

121.1 

1279. 

3  . 

13 

12 

10631.6 

34.57 

-91 . 36 

-45.3 

-53.2 

271  . 

16.7 

45.4 

0 . 0 

1168  . 

3  . 

13 

13 

10626.4 

34.55 

-91 . 22 

-45.5 

-51.0 

280  . 

16.7 

45.4 

0.0 

0  . 

0  . 

13 

14 

10625.7 

34.53 

-91.10 

-45.5 

-50.5 

293. 

19.0 

53.8 

0.0 

0  . 

0  . 

13 

15 

10620.4 

34.50 

-90.97 

-45.7 

-50.9 

298. 

17.3 

54.1 

112.3 

0  . 

0  . 

13 

16 

10613.5 

34.47 

-90.84 

-45.7 

-51.4 

305. 

15.0 

56.0 

99.2 

884. 

1  . 

13 

17 

10615.4 

34.44 

-90.69 

-45.7 

-52 . 3 

308  . 

15.0 

55.2 

114.6 

918. 

1  . 

13 

18 

10613.1 

34.41 

-90.56 

-45.9 

-53.5 

311  . 

14.7 

53.3 

119.4 

700  . 

1  . 

13 

19 

10612.0 

34.38 

-90.42 

-45.8 

-53.7 

310  . 

14.0 

55.2 

126 . 5 

1047. 

2. 
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WUSi 


13 

20 

10599.8 

34.35 

-90 . 28 

13 

21 

10618.7 

34.26 

-90 . 24 

13 

22 

10610 . 2 

34.21 

-90 . 33 

13 

23 

10628.2 

34 . 24 

-90.46 

13 

24 

10640 . 2 

34.29 

-90.57 

13 

25 

10641 . 2 

34.34 

-90.68 

13 

26 

10625.0 

34.39 

-90.80 

13 

27 

10358.4 

34.44 

-90.93 

13 

28 

10021 . 3 

34.47 

-91.04 

13 

29 

9759.9 

34.49 

-91.17 

13 

30 

9420.0 

34.51 

-91 . 29 

13 

31 

9434 . 5 

34.53 

-91.42 

13 

32 

9431.5 

34.56 

-91.53 

13 

33 

9433.3 

34.58 

-91.66 

13 

34 

9429.6 

34.61 

-91.78 

13 

35 

9429.6 

34.63 

-91 .90 

13 

36 

9418.8 

34.65 

-92.01 

13 

37 

9433.2 

34.67 

-92.13 

13 

38 

9432 . 2 

34.75 

-92 . 17 

13 

39 

9429.6 

34.79 

-92 . 06 

13 

40 

9426.6 

34.75 

-91 . 94 

13 

41 

9436.4 

34.68 

-91 .83 

13 

42 

9437.4 

34.62 

-91.74 

13 

43 

9435.4 

34 . 57 

-91 .63 

13 

44 

9427.8 

34.55 

-91 .49 

13 

45 

9405.8 

34.51 

-91 . 35 

13 

46 

9065.8 

34.48 

-91 . 25 

13 

47 

8324.0 

34.40 

-91 . 24 

13 

48 

7919.9 

34.40 

-91 . 35 

13 

49 

7930.8 

34.43 

-91.45 

13 

50 

7927.6 

34.48 

-91.54 

13 

51 

7931.4 

34.53 

-91.64 

13 

52 

7937.9 

34.57 

-91.75 

13 

53 

7940.6 

34.60 

-91.85 

33 

54 

7923.8 

34.63 

-91.97 

13 

55 

7930.5 

34.65 

-92.08 

13 

56 

7935.3 

34.67 

-92.18 

13 

57 

7936.2 

34.69 

-92.28 

13 

58 

7931.5 

34.71 

-92.39 

13 

59 

7938.8 

34.72 

-92.50 

14 

0 

7943.0 

34.75 

-92.60 

14 

1 

7935.7 

34.76 

-92.72 

14 

2 

7934.3 

34.78 

-92.83 

14 

3 

7924.3 

34.80 

-92.94 

14 

4 

7737 . 1 

34.82 

-93.04 

14 

5 

7380.5 

34.84 

-93.16 

14 

6 

7123.4 

34.86 

-93.27 

14 

7 

6537.0 

34.88 

-93.37 

14 

8 

5935.2 

34.90 

-93.47 

14 

9 

5493.3 

34.93 

-93.58 

14 

10 

5473.1 

34.95 

-93.67 

14 

11 

5469.1 

34.97 

-93.76 

14 

12 

5465.5 

34.99 

-93.86 

14 

13 

5482.1 

35.01 

-93.96 

14 

14 

5475.7 

35.03 

-94.06 

14 

15 

5475.1 

35.06 

-94.16 

14 

16 

5473.1 

35.08 

-94 . 27 

14 

17 

5479.0 

35.10 

-94.36 

14 

18 

5479.6 

35.12 

-94.46 

14 

19 

5476 . 8 

35.14 

-94.57 

14 

20 

5462.4 

35.15 

-94.69 

14 

21 

5467.1 

35.15 

-94.79 

14 

22 

5464.0 

35.15 

-94.90 

14 

23 

5479.1 

35.15 

-95.01 

14 

24 

5463.2 

35.13 

-95.11 

14 

25 

5462.6 

35.10 

-95 .22 

14 

26 

5468.6 

35.08 

-95.33 

14 

27 

5463.7 

35.06 

-95.44 

14 

28 

5471 . 3 

35.04 

-95.55 

54.1 

127 . 2 

1025 

52 . 0 

130.4 

0 

53 . 1 

123.2 

2036 

54.6 

121.9 

1940 

52.8 

0.0 

1939 

53.3 

0 . 0 

0 

54 . 3 

0 . 0 

0 

57 . 2 

130.4 

0 

58.4 

109 . 9 

0 

40.9 

106.0 

0 

68 . 2 

116.7 

68 . 3 

142 . 5 

78 . 1 

148 . 2 

86.8 

145.6 

79.7 

137.3 

I 

M. 


$ 

IS.* 

w< 

Mi 


14 

29 

5462 . 3 

35.03 

-95.67 

-6.8 

14 

30 

5467.7 

35.02 

-95.78 

-6.5 

14 

31 

5461.7 

35.01 

-95.89 

-6.5 

14 

32 

5465.0 

35.00 

-95.99 

-6.6 

14 

33 

5430.7 

34.99 

-96.11 

-6.2 

14 

34 

5208.8 

34.98 

-96 . 23 

-5 . 1 

14 

35 

4969.0 

34.98 

-96.33 

-3.5 

14 

36 

4640.4 

34.98 

-96.44 

-2.0 

14 

37 

4292.2 

35.00 

-96 . 53 

-0.2 

14 

38 

3851.2 

35.02 

-96.62 

2.9 

14 

39 

3372.4 

35.04 

-96.71 

7.2 

14 

40 

2980.6 

35.07 

-96.81 

10.1 

14 

41 

2552.2 

35.10 

-96.90 

13.8 

14 

42 

2160.9 

35 . 13 

-96.99 

15.3 

14 

43 

1728.0 

35.17 

-97.08 

16 . 5 

14 

44 

1183.4 

35.20 

-97.16 

19 . 1 

14 

45 

860 . 5 

35.23 

-97.23 

22.3 

14 

46 

862 . 2 

35.25 

-97 . 30 

22 . 3 

14 

47 

869.4 

35.27 

-97 . 38 

22 . 5 

14 

48 

870 . 0 

35.32 

-97.44 

22 .6 

14 

49 

873.6 

35.38 

-97 . 47 

23.1 

14 

50 

863.5 

35.44 

-97 . 51 

23 . 5 

14 

51 

787 . 1 

35.46 

-97.57 

24.4 

14 

52 

763.9 

35.45 

-97.62 

24.7 

14 

53 

545.4 

35.42 

-97.60 

26.7 

14 

54 

368.5 

35.40 

-97.59 

29 . 8 

14 

55 

365.4 

0.00 

0.00 

30.9 

FLIGHT 

2 

HR 

MIN 

ALT 

LAT 

LONG 

TEMP 

(M) 

(DEG  N) 

(DEG  E) 

(C) 

17 

3 

368.0 

35.40 

-97.61 

31.8 

17 

4 

368.4 

35.41 

-97.61 

31.0 

17 

5 

368.8 

35.41 

-97.60 

31.1 

17 

6 

369.8 

35.41 

-97.60 

30.7 

17 

7 

515.1 

35.39 

-97.60 

28.3 

17 

8 

1159.5 

35.35 

-97. S7 

22.3 

17 

9 

1919.1 

35.36 

-97.50 

16.1 

17 

10 

2472.3 

35.38 

-97.42 

14.1 

17 

11 

3212.1 

35.39 

-97.32 

10.2 

17 

12 

3968.3 

35.39 

-97.22 

4.1 

17 

13 

4428.1 

35.37 

-97.14 

1.6 

17 

14 

4880.0 

35.35 

-97.04 

-1.4 

17 

15 

5469 . 5 

35.34 

-96.93 

-5.3 

17 

16 

5837.9 

35.33 

-96.83 

-8.6 

17 

17 

6170.9 

35.32 

-96.74 

-11.3 

17 

18 

6616 . 3 

35.31 

-96.62 

-14.5 

17 

19 

6976.8 

35.29 

-96.50 

-17.7 

17 

20 

7001.9 

35.28 

-96.40 

-18.5 

17 

21 

6994.3 

35.27 

-96 . 28 

-18.6 

17 

22 

6989.5 

35.25 

-96.15 

-18.5 

17 

23 

6986.4 

35.24 

-96.02 

-18.4 

17 

24 

6983.7 

35.22 

-95.90 

-18.3 

17 

25 

6983.9 

35.21 

-95.78 

-18.3 

17 

26 

6989.9 

35.16 

-95.68 

-18.0 

17 

27 

6996.1 

35.08 

-95.70 

-17.7 

17 

28 

7002.3 

35.06 

-95.81 

-17.9 

17 

29 

7012 . 3 

34.99 

-95.87 

-18.1 

17 

30 

6999 . 3 

34.89 

-95.88 

-18.0 

17 

31 

7002.6 

34 . 80 

-95.90 

-17.9 

17 

32 

6999.6 

34 . 71 

-95.94 

-17.9 

17 

33 

7017.7 

34.68 

-96 . 04 

-18.2 

17 

34 

7005 . 3 

34.75 

-96.07 

-18.4 

17 

35 

7016 . 1 

34.84 

-96 . 07 

-18.5 

17 

36 

7004 . 2 

34 . 93 

-96 . 07 

-18.5 

17 

37 

7002.7 

35.01 

-96.02 

-18.6 

17 

38 

7013.3 

35.08 

-95.96 

-18.6 

17 

39 

6989.7 

35.15 

-95.89 

-18.6 

17 

40 

6978 . 9 

35.24 

-95.81 

-18.3 

-16.9 

312. 

6 . 1 

45.7 

78.9 

46  . 

1 

-17.5 

318  . 

5.7 

45.5 

80.9 

39  . 

1 

-17.3 

315. 

5.0 

46 . 1 

75.0 

0  . 

0 

-16.3 

313. 

5.7 

45.4 

81 . 3 

0  . 

0 

-17.3 

316  . 

4.7 

45.4 

87.6 

0  . 

0 

-13.6 

249. 

1.9 

48.1 

84.1 

181  . 

2 

-13.9 

281 . 

2.4 

49.6 

96.5 

169. 

2 

-14.5 

339. 

6.9 

56.9 

0.0 

277  . 

2 

-17.4 

219. 

5.8 

66 . 2 

0 . 0 

279. 

2 

-12.9 

43. 

5.1 

72.0 

0 . 0 

461 . 

2 

-8.9 

82  . 

6.0 

73.6 

121.7 

466  . 

2 

-1.7 

107  . 

7.8 

64.7 

139.6 

449  . 

2 

2 . 1 

119. 

7 . 3 

53.1 

111 .  2 

422  . 

2 

5.9 

126  . 

5 . 0 

49.4 

95.0 

580  . 

2 

9.9 

142  . 

5 . 2 

50.6 

110.7 

684. 

2 

16.1 

157  . 

6.4 

45.9 

135.7 

1087  . 

2 

18.4 

170  . 

6 . 4 

38.1 

139.8 

1371  . 

2 

18 . 2 

166  . 

7 . 0 

36.5 

117.8 

863  . 

2 

18.5 

156  . 

6.7 

45.8 

130 . 5 

670  . 

2 

19 . 5 

158. 

7 . 4 

48.5 

0 . 0 

728  . 

2 

18 . 9 

155. 

6 . 5 

47.9 

0 . 0 

0  . 

0 

18.5 

157  . 

6 . 2 

50.5 

0 . 0 

0  . 

0 

18.4 

155. 

6.0 

46.7 

0 . 0 

0  . 

0 

18.4 

155  . 

6 . 3 

41 . 5 

0.0 

0  . 

0 

18 . 5 

148  . 

7.0 

45.0 

0 . 0 

0  . 

0 

19.6 

249. 

2.9 

42.9 

0.0 

0  . 

0 

19.7 

0  . 

0.0 

44.6 

0 . 0 

0  . 

0 

TD 

WD 

WS 

03 

CO 

NO  STAT 

(C) 

(DEG) 

(M/S) 

(PPB) 

(PPB) 

(PPT) 

18.3 

0  . 

0.0 

32.8 

0.0 

0  . 

0  . 

19.3 

0  . 

0.0 

32.6 

0.0 

0  . 

0  . 

20.2 

0  . 

0.0 

32.5 

0.0 

0  . 

0  . 

20.3 

0. 

0.0 

29.0 

0.0 

0  . 

0  . 

20.2 

147. 

6.4 

22.6 

0.0 

0  . 

0  . 

18.2 

161. 

6.6 

32.9 

0.0 

0. 

0  . 

12.0 

171. 

8.6 

34.2 

0.0 

0  . 

0  . 

4.2 

164. 

5.9 

38.6 

0.0 

0  . 

0. 

-3.1 

143. 

5.3 

47.9 

0.0 

0  . 

0  . 

-4.7 

266. 

3.5 

46.7 

0.0 

0  . 

0  . 

-13.7 

341. 

9.6 

44.5 

0.0 

0  . 

0  . 

-17.3 

341 . 

8.8 

44.4 

0.0 

0  . 

0  . 

-22.0 

319. 

3.3 

33.0 

0.0 

0  . 

0  . 

-25.0 

303. 

2.9 

33.2 

92.1 

0  . 

0  . 

-27.3 

327  . 

2.8 

40 . 2 

96.0 

0  . 

0  . 

-29.3 

92  . 

4.9 

42.2 

94.5 

0  . 

0  . 

-31.2 

93. 

5.1 

37.1 

107.1 

0  . 

0  . 

-33.2 

330  . 

5.3 

43.9 

96.4 

0  . 

0  . 

-33.2 

346. 

6 . 2 

46.3 

100 . 2 

0  . 

0  . 

-29.7 

343. 

6.6 

46 . 3 

96 . 3 

0  . 

0  . 

-29.7 

340. 

6.7 

47.6 

94.8 

0  . 

0  . 

-29.9 

342  . 

6.5 

48.4 

0.0 

0  . 

0  . 

-32.4 

337. 

6.2 

44.8 

0.0 

0  . 

0  . 

-37.2 

347. 

4.9 

47.6 

0.0 

0  . 

0  . 

-39.0 

349. 

5.6 

49.7 

79.4 

93. 

1 . 

-39.4 

353. 

5.5 

50 . 0 

90 . 3 

109. 

1 . 

-39.7 

38. 

4.7 

46 . 2 

85.5 

0  . 

0  . 

-40 . 1 

16. 

4.7 

44.7 

87 . 5 

182  . 

2  . 

-40.4 

18  . 

5 . 1 

47 . 0 

91 . 5 

165. 

2  . 

-37.4 

15. 

5.8 

50 . 7 

92 . 3 

137. 

2. 

-25.6 

9  . 

5.7 

55.8 

103.0 

2727  . 

3  . 

-34.0 

18  . 

6 . 5 

46 . 2 

95.7 

1702  . 

3. 

-40 . 8 

17  . 

7.0 

44 . 0 

0 . 0 

938  . 

3  . 

-41 .4 

12  . 

6.5 

45.6 

0 . 0 

0  . 

0. 

-41.5 

8  . 

6.6 

49.1 

0 . 0 

0  . 

0  . 

-32.1 

294  . 

7.4 

49.8 

0.0 

0  . 

0  . 

-29.6 

351 . 

6.7 

50 . 1 

77 . 5 

31  . 

1  . 

-30.  J 

343. 

6.5 

49.6 

78. 0 

34  . 

1  . 
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